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ABSTRACT 

Chaff  is  the  least  expensive  and  simplest  soft  kill  measure 
against  anti-ship  missiles.  In  this  thesis  we  study  and  simulate 
the  deployment  of  chaff  in  centroid  mode  for  anti-ship  missile 
defense.  The  idea  is  to  fire  chaff  of  radar  cross  section  greater 
than  that  of  the  ship,  cause  a  shift  of  centroid  and  achieve  break 
lock  by  moving  the  ship  out  of  the  radar  resolution  cell  of  the 
missile.  After  firing  the  chaff  in  an  optimum  direction,  the  ship 
is  made  to  maneuver  in  two  steps;  first  to  make  an  earliest 
possible  break  lock  and  then  to  maximize  its  displacement  from  the 
missile  line  of  flight.  Since  the  interaction  scenario  is  time 
variant,  modelling  has  been  structured  on  a  relative  coordinate 
system.  The  algorithms  to  implement  these  objectives  are  based  on 
values  extracted  from  this  coordinate  system  centered  on  the  chaff, 
mobile  with  the  wind  and  oriented  towards  the  missile. 

The  mathematical  model  is  translated  into  a  simulation 
program.  We  consider  it  to  be  a  good  tool  for  insight  and 
understanding  and  feel  that  it  has  a  promise  of  application  both 
in  design  work  and  at  sea  in  the  passive  point  defense  problem. 


IV 


THESIS  DISCLAIMER 

The  reader  is  cautioned  that  computer  programs  developed  in 
this  research  may  not  have  been  exercised  for  all  cases  of 
interest.  While  every  effort  has  been  made,  within  the  time 
available,  to  ensure  that  the  programs  are  free  of  computational 
and  logic  errors,  they  cannot  be  considered  validated.  Any 
application  of  these  programs  without  additional  verification  is 
at  the  risk  of  the  user. 
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I .   INTRODUCTION 

A.  IMPORTANCE  OF  DEFENSE  AGAINST  ANTI-SHIP  MISSILES 

The  maritime  combat  scenario  of  today  is  filled  with  complex 
multi-dimensional  threats  for  a  man  of  war.  Conflicts  in  the  past 
quarter  century  or  so  have  shown  that  warships  have  grown  more 
vulnerable  to  aerial  threats,  in  particular  to  anti-ship  missiles. 
We  have  the  following  prominent  episodes  of  recent  times  to  endorse 
and  highlight  this  fact: 

•  An  Israeli  destroyer  EILAT  was  sunk  near  Port  Said  in  the  1967 
Arab-Israeli  war  by  Egypt  with  Styx  missile  fired  by  Soviet  built 
Komar  class  boats. 

•  The  Pakistan  Navy  destroyer  KHAIBAR  and  minesweeper  MUHAFIZ 
were  sunk  and  logistic  ship  DACCA  was  damaged  by  salvoes  of  Styx 
missiles  fired  by  Ossa  class  missile  boats,  which  are  also  Soviet 
built,  during  the  1971  Indo-Pakistan  war. 

•  Gabriel  missiles  launched  by  the  Israeli  Navy  inflicted  heavy 
damage  to  the  Syrian  Navy  resulting  in  sinking  of  a  number  of 
missile  boats  during  the  1973  Arab-Israeli  war. 

•  The  Royal  Navy  suffered  a  humiliating  blow  when  a  number  of 
its  front  line  ships  were  sunk  and  damaged  by  shore  and  air 
launched  Exocet  missiles  fired  by  the  Argentinian  Navy  during  the 
Falklands  war  in  1982  [Ref.  10]. 

•  US  ship  and  air  launched  Harpoon  missiles  destroyed  a  number 
of  Libyan  missile  boats  and  other  craft  in  March  1986  during  a 
conflict  with  Libya  in  the  Gulf  of  Sidra. 

•  In  April  1988  again,  Harpoon  and  Standard  missiles  fired  by 
the  US  Navy  caused  heavy  damage  to  Iranian  Navy  vessels  in  the 
Persian  Gulf. 

B.  APPROACHES  TO  SHIP  DEFENSE 

All  major  and  minor  maritime  forces  of  the  world  recognize  the 
vulnerability  of  surface  ships  to  this  serious  threat  today. 
However,  since  the  surface  fleets  are  foreseen  to  be  the  mainstay 
of  most  navies  for  a  long  time  to  come,  their  defense  against  anti- 
ship  missiles  has  received  tremendous  emphasis.  Larger  navies  who 


have  access  to  sophisticated  technological  and  vast  monetary 
resources  tie  their  belief  to  the  methodology  of  layered  defense 
or  defense  in  depth.  Such  a  defense  methodology  would  probably 
include  offensive  long-range  aerial  strikes  against  enemy  storage 
depots  and  platforms,  combat  air  patrols,  standoff  missiles  and  a 
variety  of  other  active  or  passive  tactics.  Even  these  elaborate 
measures  fail  to  provide  an  absolute  defense.  Smaller  navies  who 
do  not  have  such  luxuries  available  at  their  disposal  must  live 
with  lesser  and  cheaper  methods.  They  have  to  rely  on  early 
detection  of  the  enemy  and  in  most  cases  early  detection  of  the 
missile.  Such  a  restriction  imposes  the  use  of  short-range  active 
or  passive  measures.  These  measures  include  active  ECM, 
conventional  guns,  point  defense  systems  or  close  in  weapon  systems 
(CIWS)  such  as  Vulcan  Phalanx/Breda  and  chaff.  Ideally  one  could 
achieve  most  effective  results  by  having  all  of  them  working  in 
unison  in  an  intelligent  manner.  However,  availability  of  active 
ECM  and  CIWS  involves  sizable  finances  and  in  many  cases  their 
availability  is  dependent  upon  politics.  Chaff,  on  the  other  hand, 
presents  itself  to  be  a  cheap  and  simple  alternative. 

C.    MODES  OF  CHAFF  DEPLOYMENT 

Chaff  has  been  used  effectively  in  a  number  of  military 
conflicts  since  the  second  World  War.  Its  flexibility  enables  it 
to  be  used  at  various  stages  of  conflict  in  many  modes.  The 
discussion  here  will  be  restricted  to  a  confrontation  between  an 


anti-ship  missile  and  a  surface  unit  at  sea.  Chaff,  in  such  a 
situation  can  be  used  in  the  following  modes  [Ref.  13:p.  338]: 

1.  Confusion 

In  confusion  mode,  chaff  is  fired  at  long  range  (a  few 
kms)  when  the  ship  is  outside  the  range  of  enemy  search  and  target 
indication  radars.  This  tactic  provides  a  large  number  of  false 
targets  that  are  indistinguishable  from  the  real  ship. 

2 .  Distraction/Dilution 

When  the  ship  has  been  detected  by  the  enemy  and  a 
missile  attack  is  imminent,  chaff  is  fired  at  shorter  ranges  (1-2 
kms)  in  order  to  mislead  the  missile  radar  when  it  starts 
transmitting.  The  chaff  cloud's  radar  cross  section  need  only  be 
comparable  to  the  ship  in  this  mode. 

3 .  Dump 

This  mode  is  also  called  range  gate  pull  off/stealing. 
Once  the  missile  is  locked  on,  chaff  is  fired  at  short  ranges  (less 
than  1  km) .  In  this  mode  the  ship's  ECM  equipment  is  also  used  in 
conjunction  with  the  chaff.  The  ship's  jammer  misleads  the  missile 
on  to  the  chaff  and  when  jamming  ceases,  the  missile  homes  onto  the 
chaff  decoy.  In  other  words,  break  lock  is  achieved  with  the  help 
of  the  jammer. 

4 .  Cent ro id 

In  this  mode,  chaff  is  deployed  at  very  close  ranges. 
The  exact  figure  depends  on  the  philosophy  of  actual  deployment. 
Some  experts  suggest  firing  chaff  at  a  range  of  25  to  50  m  so  that 
its  echo  merges  with  that  of  the  ship  and  the  missile  homes  onto 


the  missile  is  lured  away  avoiding  a  direct  hit  thus  minimizing 
damage.  This  is  a  last  ditch  method  as  break  lock  is  not  the 
primary  goal. 

D.  RECOMMENDED  SOLUTION 

Although  reaction  times  are  shrinking  with  sophisticated 
missiles  entering  the  battle  arena,  intelligent  deployment  of  chaff 
has  been  known  to  present  an  effective  counter.  Ideally,  chaff 
could  be  most  effective  if  integrated  with  other  forms  of  active 
and  passive  anti-ship  missile  defenses.  However,  if  monetary 
constraints  preclude  such  a  luxury,  we  consider  that  chaff  provides 
a  good  alternative  to  thwart  most  missile  threats  of  today.  It  can 
be  argued  that  recent  advances  in  technology  have  enabled 
production  of  smart  and  sophisticated  missile  radars  which  have  the 
ability  to  reject  chaff  clutter  by  sensing  the  difference  of  its 
parameters  from  those  of  the  ship,  thereby  rendering  it 
ineffective.  Nevertheless,  there  is  always  another  counter  to  a 
counter.  Chaff  systems  are  also  getting  smarter  with  added 
features.  These  features  include  chaff's  maneuvering  capability 
with  the  help  of  a  self  propulsion  system  inside  the  chaff  and 
expendable  ECM  devices  with  integrated  processing  facility.  Active 
decoy  developments  such  as  Nulka  and  Winnin  are  examples  of  such 
chaffs . 

E.  AIM  OF  THE  THESIS  AND  RELATED  WORK 

The  purpose  of  this  work  is  to  simulate  centroid  mode  of  chaff 
deployment,  fire  the  chaff  at  close  range  (100  to  400  m) ,  achieve 


break  lock  and  maximize  miss  distance  by  intelligent  ship  maneuvers 
assuming  a  supporting  jammer  is  not  available.  Work  has  also  been 
done  in  this  field  by  David  Nitzan  and  Chris  Nygard  [Refs.  22  &  23] 
The  mathematical  modelling  and  simulation  done  in  this  thesis  has 
been  developed  as  an  extension  to  their  work. 


II.   KINEMATICS  OF  THE  PROBLEM 

A.    MISSILE  DETECTION 

One  of  the  most  formidable  and  decisive  problems  in  the  field 
of  anti-missile  defense  is  early  detection  of  the  missile.  Most  of 
today's  anti-ship  missiles  travel  at  speeds  close  to  MACH  1  and  fly 
very  low  to  the  surface  of  the  sea.  This  sea  skimming  capability 
coupled  with  the  state-of-the-art  aerodynamic  characteristics, 
makes  the  task  of  missile  detection  even  more  difficult  for  the 
Operations  Room  of  the  ship.  For  vessels  operating  without  any  air 
support,  reliance  has  to  be  placed  on  sensors  available  on  board. 
The  two  key  sensors  of  help  in  this  situation  are  radar  and  ESM 
equipment.  ESM  can  be  coupled  with  radar  to  confirm  detection  made 
difficult  by  background  clutter  at  low  altitudes.  However,  modern 
missiles  commence  radiating  for  their  terminal  homing  phase  quite 
late  in  their  flight,  thus  limiting  the  ESM  capability  to  ranges 
less  than  the  radiation  commencement  zone. 

The  primary  sensor,  as  a  matter  of  principle,  should  be  radar. 
Nonetheless,  because  of  low  flying  profiles  of  the  missile  and  non- 
availability of  state-of-the-art  radars  to  smaller  navies,  ESM 
takes  over  the  primary  sensory  role.  This  is  especially  true  for 
the  sea  skimmers.  The  range  at  which  the  missile  commences  its 
terminal  phase  has  a  very  large  variability,  depending  on  the  type 
and  generation  of  the  missile  technology.  Based  on  available 
unclassified  data,  this  range  varies  from  around  3  miles  to  20 
miles.  With  an  average  speed  of  600  knots,  at  which  most  missiles 


travel,  we  get  a  reaction  time  ranging  from  18  seconds  to  about  2 
minutes.  It  is  quite  obvious  that  reaction  time  is  the  key  factor 
in  today's  warfare.  With  the  missiles  becoming  supersonic  and  their 
terminal  phases  rapidly  shortening,  reaction  time  is  increasingly 
diminishing.  Table  1  shows  reaction  times  for  different  initial 
detection  ranges  and  is  a  useful  indicator  to  emphasize  this  point. 
Assuming  the  missile  speed  to  be  600  knots,  the  following  data  has 
been  tabulated: 

TABLE  1 
MISSILE  RANGE  AND  REACTION  TIME 


S.No 

Missile 

Range (yds) 

Reaction 

Time (sees) 

1. 

20,000 

60 

2. 

15,000 

45 

3. 

14,000 

42 

4. 

13,000 

39 

5. 

12,000 

36 

6. 

11,000 

33 

7. 

10,000 

30 

8. 

9000 

27 

9. 

8000 

24 

10. 

5000 

15 

11. 

3000 

9 

12. 

1000 

TOO 

LATE 

B.    CHAFF  DEPLOYMENT  IN  CENTROID  MODE 

There  are  two  philosophies  associated  with  centroid  or 
seduction  mode  tactic;  the  first  is  to  fire  the  chaff  at  very  close 
ranges,  i.e.,  25  to  50  meters  [Ref.  5],  within  the  radar  resolution 
cell  of  the  missile.  The  missile  takes  the  two  echoes  to  be  one  and 


locks  on  to  the  centroid  of  the  two.  The  ship  then  moves  away 
briskly  from  the  chaff,  depending  on  the  wind  direction  leaving  the 
missile  to  home  on  the  chaff.  Depending  on  all  the  variables 
involved,  the  ship  can  either  break  lock  if  it  is  out  of  the 
resolution  cell  in  time  or  else  can  maximize  displacement  from  the 
missile  impact  point.  In  this  way,  the  ship  can  either  escape 
damage  altogether  or  at  least,  minimize  it.  Real  time  data  from 
relevant  conflicts  is  not  available  to  seek  a  recommended  figure 
for  safe  distance  of  the  ship  from  the  point  of  impact.  However, 
at  least  50  to  100  meters  is  considered  to  be  a  safe  distance  for 
minimal  damage  to  the  ship. 

The  second  philosophy  behind  utilization  of  chaff  in  centroid 
mode  is  to  fire  the  chaff  within  the  resolution  cell  on  a  bearing 
related  to  the  wind  direction  and  speed.  Chaff  is  placed  as  far  as 
possible  from  the  ship  towards  the  boundaries  of  the  resolution 
cell  [Ref.  13],  depending  on  the  initial  detection  range  of  the 
missile,  beam  width  of  the  missile  radar  and  the  wind  parameters. 
Once  the  chaff  has  bloomed,  the  missile  amends  its  homing  on 
direction  to  the  centroid.  In  the  meantime,  concurrent  to  chaff 
launch,  the  ship  maneuvers  its  way  out  of  the  resolution  cell  to 
achieve  a  break  lock  and  maximize  distance  from  the  missile  impact 
point.  Analysis  of  this  philosophy  is  the  aim  of  this  thesis. 

C.    RADAR  RESOLUTION  CELL 

The  resolution  cell  is  the  central  factor  of  the  anti-ship 
missile  problem.  Basically,  the  radar  has  a  certain  range  and 


bearing  resolution  based  on  its  effective  pulse  width  and  beam 
width  [Ref.  3:p.  332].  The  area  outside  these  bounds  is  the  limit 
of  the  radar  to  distinguish  one  target  from  another.  If  there  are 
two  targets  and  the  radar  is  locked  onto  one  of  them,  the  second 
should  be  outside  the  boundary  to  be  discriminated  as  a  separate 
target.  If  the  second  target  also  lies  within  this  area,  the  radar 
will  only  see  them  as  one  target  and  align  itself  to  their 
centroid,  both  radially  and  azimuthally.  Figure  1  shows  what  a 
radar  resolution  cell  looks  like. 


Beamwidth     O 


Missile 


ACTUAL 


Range  Gate 


Angle  Gate 


Figure  1.  Radar  Resolution  Cell  [Ref  23] 

The  limits  in  both  angle  and  range  are  defined  by  the  half 
power  or  3  dB  values.  Most  missile  seekers  use  this  geometry  for 
processing  purposes.  From  the  missile's  point  of  view,  only  the  RF 
energy  inside  this  cell  is  considered  to  be  the  useful  backscatter. 
Any  return  from  the  target  outside  this  area  is  ignored  to  ensure 
sharp  resolution.  The  cell  actually  assumes  an  annular  shape,  but 
at  ranges  we  are  interested  in,  it  can  be  approximated  to  be  a 


rectangle  without  incurring  any  sizable  inaccuracies  as  shown  in 
Figure  2 : 


Beamwidth 


Missile 


APPROXIMATION 


Figure  2.  Approximation  to  the  Resolution  Cell  [Ref.  23] 

The  bounds  of  the  cell  are  also  called  gates;  the  one 
perpendicular  to  missile  direction  is  the  angle  gate  and  the  one 
parallel  to  it  is  the  range  gate  as  shown  in  Figure  3. 
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Figure  3.  Range  and  Angle  Gate 
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The  center  of  the  cell  is  the  centroid  for  the  missile  radar 
to  lock  on.  As  the  bearing  of  the  missile  changes  due  to  movement 
of  the  ship  or  centroid,  the  cell  rotates  to  reorient  itself  in 
order  to  maintain  orthogonality  with  the  missile  line  of  flight. 
Similarly,  due  to  variation  of  target  range  from  the  missile  as  it 
closes  in,  the  dimensions  of  the  cell  vary  in  the  azimuthal  plane. 
The  radar  resolution  cell,  therefore,  forms  an  orthogonal 
coordinate  system  which  is  varying  with  time  and  constantly 
readjusting  its  bounds.  This  motion,  as  mentioned,  depends  on  the 
movement  of  the  centroid  and  the  longitudinal  axis  of  the  missile. 

D.    EFFECT  OF  CHAFF  ON  RADAR  RESOLUTION  CELL 

The  whole  idea  of  using  chaff  in  a  favorable  manner  is  to 
shift  the  centroid  of  the  resolution  cell  away  from  the  ship,  walk 
it  off  to  the  chaff  and  move  out  of  the  cell  to  break  lock  thereby 
seducing  the  missile  onto  the  chaff  cloud.  The  chaff,  of  course, 
has  to  be  fired  to  bloom  inside  the  resolution  cell.  If  the  radar 
cross  section  of  the  chaff  is  comparable  to  that  of  the  ship,  the 
missile  is  expected  to  shift  its  lock  to  the  centroid  of  the  two 
radar  cross  sections  which  is  approximately  the  mid-point  on  the 
line  joining  the  ship  and  the  chaff  [Ref.  13:p.  347]. 

Usually,  most  captains  would  like  to  fire  chaff  with  a  radar 
cross  section  considerably  larger  than  that  presented  by  the  ship. 
The  shift  of  the  centroid  would  then  be  more  favorable  to  the  ship 
and  closer  to  the  chaff  cloud.  In  any  case,  as  soon  as  the  chaff 
blooms,  it  starts  moving  with  the  wind.  The  aim  is  to  fire  the 
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chaff  on  a  bearing  that  optimally  assists  the  ship  for  the  wind 
vector  existing  during  interaction.  The  job  of  the  ship  then  is  to 
move  out  of  the  resolution  cell,  in  a  direction  that  will  take  it 
away  from  the  chaff. 

The  objective  of  moving  the  ship  out  of  the  cell  is  to  reduce 
and  then  eliminate  its  radar  cross  section  from  the  bounds  of  the 
cell.  As  the  ship  starts  moving  away  from  the  chaff,  its  distance 
from  the  centroid  also  increases.  As  it  crosses  the  cell  boundary, 
either  through  its  range  or  angle  gate,  its  radar  cross  section 
inside  the  cell  decreases.  The  effect  of  this  development  is  two- 
fold: first  the  ship  has  moved  away  from  the  centroid;  second,  due 
to  its  decreased  radar  cross  section  inside  the  cell,  the  centroid 
moves  further  towards  the  chaff.  This  effect  multiplies  as  more  of 
the  ship  exits  the  cell  and  hastens  the  shift  of  centroid  to  the 
chaff.  As  the  ship  leaves  the  cell  completely,  break  lock  is 
achieved.  The  missile  now  locks  onto  the  chaff  and  the  ship 
endeavours  to  maximize  its  distance  from  the  new  centroid,  i.e., 
the  chaff. 

The  kinematics  are,  however,  not  as  simple  as  they  appear  to 
be.  The  most  crucial  factor  at  this  stage  is  the  wind  speed  and 
direction  in  relation  to  both  the  missile  line  of  flight  and  the 
ship's  initial  course.  The  wind  could  be  a  helping  hand  or  prove 
unfriendly  depending  on  what  it  is  doing  during  the  operation.  If 
it  is  strong  and  moving  parallel  to  the  missile  line  of  flight, 
i.e.,  perpendicular  to  angle  gate,  the  chaff  will  spend  very  little 
time  in  the  cell.  This  time  might  not  be  enough  to  enable  the  ship 
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to  exit  the  cell  before  the  chaff  does  so.  As  a  result,  the  missile 
locks  on  to  the  ship  again. 

As  an  exception,  the  ship  and  the  chaff  might  decide  to  leave 
the  boundaries  of  the  cell  simultaneously.  In  such  a  situation,  it 
is  hard  to  tell  what  the  missile  will  lock  on  to.  There  appears  to 
be  an  equal  probability  of  the  missile  either  locking  on  to  the 
ship  or  the  decoy.  However,  it  first  depends  on  the  decision-making 
algorithm  of  the  processor  of  the  missile's  radar.  It  would  also 
depend  on  the  relative  location  of  the  chaff  and  the  ship  with 
respect  to  the  missile.  If  the  processor  unit  conducts  a  far-to- 
near  processing  pattern,  the  nearer  of  the  two  will  be  selected. 
Of  course,  the  reverse  would  be  true  for  a  near-to-far  pattern. 

In  summary,  the  presence  of  chaff  redefines  the  boundaries  of 
the  radar  resolution  cell,  reorients  it,  displaces  it  and  aids  the 
ship  to  get  out  of  troubled  waters.  The  crucial  factor  is  a 
favorable  wind  vector  at  the  scene  of  the  battle. 
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III.   MATHEMATICAL  MODELLING 

In  order  to  formulate  a  simulation  program,  it  is  necessary 
to  construct  a  mathematical  model  of  the  situation  in  which  the 
chaff  is  to  be  deployed.  There  are  a  number  of  parameters  that  are 
made  available  to  us  by  the  ship,  missile,  wind  and  the  chaff.  Some 
of  these  parameters  are  constant  and  others  are  variable.  The 
variables  are  classified  into  two  categories;  those  that  are 
presented  to  us  by  a  particular  situation  or  the  input  variables, 
and  those  that  we  have  the  ability  to  control  and  manipulate.  It 
is  the  latter  type  of  variables  that  make  it  possible  for  us  to 
successfully  use  the  chaff.  A  list  of  input  and  control  variables 
are  given  in  Tables  2  and  3.  The  abbreviations  for  all  the 
variables  in  this  chapter  are  chosen  to  correspond  to  those  used 
in  the  simulation  program. 

A.    SHIP'S  EQUATIONS 

As  summarized  in  the  tables,  the  ship's  parameters  consist  of 
its  length,  breadth,  course,  speed,  acceleration,  turn  rate  and 
radar  cross  section. 

1.    Initial  Conditions 

The  ship  is  assumed  to  be  at  the  center  of  the  reference 
coordinate  system  for  the  initial  set  up  at  the  time  the  missile 
is  detected.  This  is  actually  the  screen  center.  The  initial 
velocity  vector  of  the  ship  is  determined  by  its  initial  course  and 
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speed. The  following  equations  describe  the  ship  at  the  starting 

point: 

XS(h)  =  XCENT  (3.1) 

YS(h)  =  YCENT  (3.2) 

XSD(h)  =  VS(h)*Sin{ (CSE(h) **/180}  (3.3) 

YSD(h)  =  VS(h)*Cos{ (CSE(h) **/180)  (3.4) 

TABLE  2 
LIST  OF  INPUT  VARIABLES 


S.No. 

VARIABLE 

ABBREVIATION 

MISSILE 

1. 

Initial  Bearing 

BRGM(l) 

2. 

Initial  Range 

RANGE 

3. 

Missile  Speed 

VMM 

4. 

Radar  Beamwidth 

e 

5. 

Range  Gate 

RGATE 

6. 

Missile  Acceleration 

AM 

7. 

Proportional  Navigation  Constant 

N 

SHIP 

8. 

Course 

CSE(l) 

9. 

Speed 

VS 

10. 

Length 

SL 

11. 

Breadth 

SB 

12. 

Turn  Rate 

TR 

13. 

Acceleration (Velocity  Rate) 

VR 

14. 

Max  Radar  Cross  Section 

RCSSMAX 

15. 

Min  Radar  Cross  Section 

RCSSMIN 

CHAFF 

16. 

Radius 

CRAD 

17. 

Radar  Cross  Section 

RCSC 

WIND 

18. 

Wind  Direction 

GW 

19. 

Wind  Speed 

VW 

MISCELLANEOUS 

20. 

Screen  Scale 

SC 
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TABLE  3 
CONTROL  VARIABLES 


S.No. 

VARIABLE 

ABBREVIATION 

SHIP 

1. 

First  Leg  of  New  Course 

CSE2 

2. 

Second  Leg  of  New  Course 

CSE3 

3. 

Max  Speed 

VSMAX 

CHAFF 

4. 

Bearing  to  launch  chaff 

BRGC 

where 

XS(h)  and  YS(h)  are  coordinates  of  ship  displacement 
XSD(h)  and  YSD(h)  are  coordinates  of  ship  velocity 
XCENT  and  YCENT  are  coordinates  of  screen  center 
CSE(h)  is  the  ship's  course 
2.    Ship's  Movement 

Having  detected  the  missile,  the  ship  after  firing  the 

chaff,  accelerates  and  maneuvers  out  of  the  radar  resolution  cell. 

In  order  to  update   its  position,   velocity  and  acceleration 

components,  the  following  equations  are  used: 

XS(h+l)  =  XS(h)  +  XSD(h)  +  0.5  XSDD(h)*h2  (3.5) 

YS(h+l)  =  YS(h)  +  YSD(h)  +  0.5  YSDD(h)*h2  (3.6) 

where 

XSD(h)  =  VS*Sin{CSE(h)*7r/180}  (3.7) 

YSD(h)  =  VS*Cos{CSE(h)*7r/180}  (3.8) 

The  velocity  components  of  the  ship  are  shown  in  Figure  4. 
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Figure  4.  Ship's  Velocity , Components 

Also 

VS(h+l)  =  VS(h)  +  VR*h  (3.9) 

XSDD(h)  =  -VS(h)  (TR*jr/180)  *  U  *  Sin(CSE  (h)  **/180 } 

+  VR*Sin{CSE(h)**/180)  (3.10) 

YSDD(h)    =    -VS(h)  {TR*?r/180)    *    U    *    Cos  {  CSE  (h)  *tt/180  ) 

+    VR(h) *Cos{TR*x/180)  (3.11) 

We  have  supposed  that  VS  can  attain  a  maximum  value  of  25 
knots.  The  terms  used  in  the  above  equations  are  defined  as 
follows: 

XSDD(h)  and  YSDD(h)  =  Instantaneous  coordinates  of  ship's 

acceleration 

XSD(h)  and  YSD(h)  =  Instantaneous  coordinates  of  ship's 

velocity 

VS  =  Velocity  of  ship 

VR  =  Velocity  rate  or  acceleration  of  the  ship 

U  =  Turn  direction  (discussed  later) 

h  =  Time  step  increment 
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Updating  the  bearing  of  the  missile  from  the  ship  is  carried 
out  using  the  following  equations: 


BRGM(h)  =  180/*  *  tan 


-1 


XM(h+l)  -  YSfh-H) 
YM(h+lj  -  Ys(h+lj 


(3.12) 


where  XM(h+l)  and  YM(h+l)  are  updated  coordinates  of  the  missile 
position. 

3 .    Radar  Cross  Section 

The  radar  cross  section  of  the  ship  is  a  variable 
dependent  upon  the  aspect  presented  by  the  ship  to  the  missile 
radar.  It  is  assumed  that  the  missile  presents  maximum  radar  cross 
section  when  viewed  from  right  ahead,  right  astern  and  the  beams. 
As  the  aspect  changes,  the  radar  cross  section  decreases  to  a 
minimum  at  the  bows  and  the  quarters.  In  the  absence  of  real  time 
data,  we  consider  this  assumption  to  be  a  reasonable  approximation. 
It  is  modelled  with  the  help  of  the  following  equation: 

RCSS  =  RCSSMIN+ |RCSSMAX-RCSSMIN*Cos [2 (CSE (1) -BRGM(l)  }*A/180]  | 

(3.13) 
where 

RCSS  =  Radar  cross  section  of  the  ship 

RCSSMAX  =  Maximum  radar  cross  section  of  the  ship 

RCSSMIN  =  Minimum  radar  cross  section  of  the  ship 

CSE(l)  =  Initial  course  of  the  ship 

BRGM(l)  =  Initial  bearing  of  the  missile 
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Pictorially,  the  equation  produces  a  graph  of  the  ship's  radar 
cross  section  vs  the  aspect  angle  as  shown  in  Figure  5. 
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Figure  5.  Radar  Cross  Section  Model 
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4.    The  Ship's  Model 

The  ship's  geometry  is  approximated  by  a  polygon  centered 
on  (XS,YS)  using  its  length  and  breadth.  The  orientation  of  the 
ship,  pivoted  on  its  center  is  shown  with  the  help  of  its  course 
with  respect  to  true  North  as  shown  in  the  Figure  6. 

The  equations  to  find  its  vertices  are  given  as  under: 
XSHIP(l)  =  XS(h)  +  SL/2  *  Sin  (CSE(h)}  (3.14) 

YSHIP(l)  =  YS(h)  +  SL/2  *  Cos  (CSE(h)}  (3.15) 

SHIPS  MODEL 


North 


CSE 


Xship(5),   Yship(5) 


Xship(4),   Yship(4) 


Xship(1),   Yship(1) 


Xship(2),   Yship(2) 


Xship(3),   Yship(3) 


Figure  6.  The  Ship's  Model  [Ref.  23] 
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XSHIP(2)    =    XS(h)    +    SL/4    *    Sin    <CSE(h)}    +    SB/2    *    Cos{CSE(h)} 

(3.16) 
YSHIP(2)    =   YS(h)    +    SL/4    *    Cos    (CSE(h)}    -   SB/2    *    Sin    {CSE(h)} 

(3.17) 
XSHIP(3)    =   XS(h)    -   SL/2    *    Sin    {CSE(h)}    +   SB/2    *    Cos    {CSE(h)} 

(3.18) 
YSHIP(3)    =   YS(h)    -   SL/2    *    Cos    (CSE(h)}    -   SB/2    *    Sin    {CSE(h)} 

(3.19) 
XSHIP(4)    =    XS(h)     -    SL/2    *    Sin    (CSE(h)}    -    SB/2    *    Cos    {CSE(h)) 

(3.20) 
YSHIP(4)    =   YS(h)    -   SL/2    *    Cos    (CSE(h)}    +   SB/2    *    Sin    {CSE(h)} 

(3.21) 
XSHIP(5)    =    XS(h)    +    SL/4    *    Sin    (CSE(h)}    -    SB/2    *    Cos    (CSE(h)} 

(3.22) 
YSHIP(5)    =    YS(h)    +    SL/4    *    Cos    (CSE(h)}    +    SB/2    *    Sin    (CSE(h)} 

(3.23) 

B.    CHAFF  EQUATIONS 

The  following  parameters  are  related  to  the  chaff: 

•  Launch  bearing 

•  Distance  from  the  ship 

•  Radius  of  the  chaff 

•  Radar  cross  section  of  the  chaff 

•  Bloom  time 

•  Velocity  vector  of  chaff  after  bloom 

1.    Initial  Conditions 

Since  the  wind  is  the  fundamental  issue  here,  it  is 
appropriate  to  see  its  relationship  with  the  chaff.  Wind  direction 
is,  by  convention,  taken  to  be  the  direction  it  is  blowing  from, 
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the  direction  in  which  the  chaff  moves  is  its  reciprocal.  This  is 
shown  in  Figure  7.  The  unsealed  velocity  of  the  wind  is  VW1  and  it 
is  also  the  velocity  of  the  chaff. 


Velocity    Vector 


Wind 
Dir. 


GC=GW+180 


VW1=VC 

Figure  7.  The  Wind  Velocity  Vector 

where 

VW1  =  Wind  speed  (Unsealed) 

GW  =  Wind  direction 

Velocity  vector  of  the  chaff  is  given  by  the  following 
equation: 

VC  =  VW 

GC  =  GW  +  180  (3.24) 

Let  us  now  consider  the  initial  positioning  of  the  chaff.  In 
order  to  offset  the  centroid  as  far  away  as  possible,  we  fire  the 
chaff  so  that  its  center  lies  on  the  boundary  of  the  angle  gate. 
The  angle  gate  is,  however,  shrinking  as  the  missile  closes  in. 
The  initial  placement  of  the  chaff,  therefore,  is  brought  closer 
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to  the  cell  center  laterally  by  one  chaff  radius  in  order  to  make 
it  bloom  inside  the  future  dimensions  of  the  cell.  The  calculation 
of  this  distance  is  done  using  the  following  equation: 

CDIST  =  R(l)  *  TAN(THETA  *  PI/180)  -  CRAD  (3.25) 

The  initial  position  of  the  chaff  round  is  shown  in  Figure  8. 


North 


Launching  Vector 


(XC,  YC) 
BRGC  ^ 


Ship 


CDIST 


CDIST  Cos  BRGC 


CDIST  Sin  BRGC 
(XCENT,  YCENT)  =  (XS,  YS) 


Figure  8.  Initial  Position  of  Chaff  [Ref.  23] 

The  parameters  in  the  figure  are  as  follows: 

BRGC   =  Chaff  launch  bearing 

CDIST  =  Distance  of  chaff  from  the  ship 

We  can  see  here  that  the  coordinates  of  the  chaff  are: 

XC(h)  =  XCENT  +  CDIST  *  Sin(BRGC  *  jt/180)  (3.26) 

YC(h)  =  YCENT  +  CDIST  *  Cos(BRGC  *  tt/180)  (3.27) 

2.    Chaff  Movement 

The  ship  is  initially  in  the  screen  center.  As  the  chaff 

starts  blowing  with  the  wind,  its  coordinates  are  updated  by  the 

following  equations: 
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XC(h+l)  =  XC(h)  +  VC  *  Sin  (GC  *  */180)  (3.28) 

YC(h+l)  =  YC(h)  +  VC  *  Cos  (GC  *  */180)  (3.29) 

3.  Bloom  Time  Delay 

In  order  to  Maximize  the  factor  of  realism,  we  have  taken 
into  account  the  time  delay  between  the  firing  of  the  chaff  and  its 
full  bloom.  For  simplicity,  it  has  been  assumed  that  the  time  taken 
to  launch  the  chaff  is  included  in  the  bloom  time  delay.  Although 
various  types  of  chaff  have  different  time  delays  in  this  context, 
we  consider  a  bloom  time  delay  of  5  seconds  as  an  appropriate 
average  figure  [Ref.  13:p.  342].  It  also  includes  the  time  it  takes 
to  launch  the  chaff  at  its  intended  initial  position.  The 
implmentation  of  this  concept  is  achieved  by  keeping  the  chaff  on 
the  ship  for  the  first  5  iterations.  On  the  sixth  second,  the  chaff 
cloud  comes  into  existence  at  the  desired  location.  Thereafter  it 
is  at  the  mercy  of  the  wind  vector  and  starts  drifting.  Since  the 
ship  and  the  resolution  cell  are  mobile,  there  is  a  possibility 
that  the  chaff  may  bloom  outside  the  cell  as  a  result  of  bloom  time 
delay.  To  avoid  such  an  eventuality,  the  first  position  of  the 
bloomed  chaff  is  offset  in  accordance  with  the  ship's  movement. 

4.  The  Chaff  Model 

The  chaff  cloud  is  approximated  by  a  circle  but  due  to 
convenience  of  mathematical  calculations,  it  is  drawn  as  a  hexagon 
using  graphics.  The  hexagon  is  centered  on  the  chaff  center  (XC,YC) 
and  is  shown  in  Figure  9.  Coordinates  for  the  corners  are  worked 
out  with  the  help  of  the  following  equations: 
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XCF(l)    =    XC(h)    +    CRAD 
YCF(l)    =    YC(h) 

XCF(3),  YCF(3) 


XC.YC 


XCF(2),  YCF(2) 


XCF(4),  YCF(4)  fa    CRAp 


XCF(1),  YCF(1) 


XCF(5),  YCF(5) 


XCF(6),  YCF(6) 


XCF(2) 
YCF(2) 
XCF(3) 
YCF(3) 
XCF(4) 
YCF(4) 
XCF(5) 
YCF(5) 
XCF(6) 
YCF(6) 


Figure  9.  Chaff  Cloud  Model  [Ref.  23] 

XC(h)  +  CRAD  *  Cos (60°) 

YC(h)  +  CRAD  *  Sin(60°) 

XC(h)  -  CRAD  *  Cos(60°) 

YC(h)  +  CRAD  *  Sin(60°) 

XC(h)  -  CRAD 
YC(h) 

XC(h)  -  CRAD  *  Cos(60°) 

YC(h)  -  CRAD  *  Sin(60°) 

XC(h)  +  CRAD  *  Cos(60°) 

YC(h)  -  CRAD  *  Sin(60°) 


(3.30) 
(3.31) 


(3.32) 
(3.33) 
(3.34) 
(3.35) 
(3.36) 
(3.37) 
(3.38) 
(3.39) 
(3.40) 
(3.41) 


C.    EQUIVALENT  TARGET  EQUATIONS 

Before  the  chaff  is  fired,  the  missile  radar  is  locks  on  to 
the  ship  and  is  tracking  it.  When  the  chaff  is  fired  and  it  has 
bloomed  inside  the  radar  resolution  cell  of  the  missile,  the  radar 
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senses  a  new  situation.  The  weighted  center  of  the  two,  i.e.,  the 
ship  and  the  chaff,  is  now  seen  by  the  missile  as  the  centroid  of 
the  target  and  it  starts  tracking  it.  This  centroid,  also  called 
the  equivalent  target,  is  determined  by  the  magnitude  of  the  radar 
cross  sections  of  the  chaff  and  the  ship  as  seen  by  the  missile. 
It  is  calculated  by  the  following  equations: 

vira'/ii.i.il     XSjhJ  *RCSS(h)  +  XC(h)  *  RCSC(h)  ,*    .,> 

XET(h+l)  -  RCSSfhj  +  RcSC(h) (3.42) 

VFT,h+n     YS(h)  *RCSS(h)  -t-  YC(h)  *  RCSC(h)  ,^    Ao* 

YET(h+l)  -       RCSSfhj  +  RcSC(h)  (3.43) 

The  velocity  with  which  this  equivalent  target  moves  is  calculated 
using  the  following  equations: 


VET  =  [(XET(h)  -  XET(h-l)}2  +  (YET(h)  -  YET (h-1) }*]*    (3.44) 


TFT  -  -t-ar.-1  .     YET(h)  -  YET(h-l) 
GET  -  tan   «  XETjhi  -  XET(h-lj 

This  is  shown  in  Figure  10. 

D.    MISSILE  EQUATIONS 

The  missile  parameters  consist  of  the  following 

Range  and  bearing  from  the  ship 

Range  resolution 

Missile  speed 

Half  power  beamwidth  of  the  antenna 

Proportional  navigation  constant 

Acceleration  of  the  missile 

Turning  rate 
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(3.45) 


{XET(h),YET(h)} 


{XET(h-1)},YET(h-1)} 

Figure  10.  Equivalent  Target  Velocity 
1.    Initial  Conditions 

The  first  parameters  related  to  the  missile  that  are 
observed  by  the  ship  are  its  initial  detection  range  and  bearing 
which  are  shown  in  Figure  11.  The  related  equations  are  as  follows: 
XM(i)  =  XCENT  +  R(i)*  Sin{BRGM(i)  *  tt/180}  (3.46) 

YM(i)  =  YCENT  +  R(i)*  Cos{BRGM(i)  *  tt/180)  (3.47) 

The  direction  flight  of  the  missile  is  the  reciprocal  of  its 
bearing  from  the  ship  i.e., 

GM  =  BRGM(i)  +18  0  (3.48) 

The  initial  velocity  components  of  the  missile  in  the  X  and  Y 
directions  are  shown  in  Figure  12.  The  eguations  for  these 
components  are  as  follows: 

XMD(h)  =  VM  *  Sin(GM  *  */180)  (3.49) 

YMD(h)  =  VM  *  Cos(GM  *  jr/180)  (3.50) 
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(XCENT,  YCENT) 


Missile 

(XM,  YM) 


R(h)  Cos  BRGM(h) 


R(h)  Sin  BRGM(h) 


Figure  11.  Missile  Range  and  Bearing 


YMD(h) 


GM(h) 

(XM,  YM) 


Figure  12 .  Missile  Velocity  Components 
2.    Missile  Movement 

The  coordinates  of  the  missile  position  as  it  moves  are 
updated  as  follows: 
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XM(h+l)  -  XM(h)  +  XMD(h)  *  h  +  0.5  *  XMDD(h)  *  h2      (3.51) 
YM(h+l)  -  YM(h)  +  YMD(h)  *  h  +  0.5  *  YMDD(h)  *  h2       (3.52) 
XMD  and  YMD  are  extracted  from  the  equations  in  the  previous 
section  whereas,  XMDD  and  YMDD,  which  represent  the  components  of 
the  missile  acceleration  are  shown  in  Figure  13  and  given  by  the 
the  following  equations: 

XMDD  =  AM  *  Sin(GM)  (3.53) 

YMDD  =  AM  *  Cos(GM)  (3.54) 


(XCENT,  YCENT) 


Missile 


(XM,  YM) 


YMDD 


XMDD 


Figure  13:  Missile  Acceleration  Components 

The  missile  direction  angle  is  updated  using  the  following 
equation: 


PM  -  fan"1    XM(h+l)  -  XM(h) 
GM  "  tan   I  YM(h+lj  -  YM(hj 


(3.55) 
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This  is  shown  in  Figure  14.  The  X  and  Y  components  of  the  missile 
position  vector  are  updated  iteratively. 


(Xm(h+1),  Ym(h+1)) 


AV 


(XM(h),  YM(h)) 


YM(h+1)  -  YM(h) 


XM(h+1)  -  XM(h) 


Figure  14.  Updating  Missile  Direction 
3.    Target  Tracking  by  the  Missile 

In  order  to  maintain  track  of  the  target  the  missile  uses 
the  principle  of  proportional  navigation  constant  i.e  it  maintains 
a  constant  line  of  sight  with  respect  to  a  reference.  This  angle 
is  shown  in  Figure  15, 
where 


+-»n'1    YET(h)  -  YM(h) 
°    ~    tan   1  XETJhj  -  XM(hj 


(3.56) 


In  order  to  maintain  a  constant  line  of  sight  angle,  the 
missile  guidance  system  senses  the  rate  of  change  in  o.  This  is 
given  by  the  following  equation: 

d/dt(a)  =  SIGD  =  VET*Sin{  (GET-a  )  ^-  (VM*Sin(  (GM-a)  )      (3>57) 
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(YET-YM) 


REFERENCE 


(XET-XM) 


Figure  15.  Target-Missile  Line  of  Sight  [Ref.  23] 

where 

RMT(h)  =  Range  of  the  missile  from  equivalent  target 

In  order  to  cater  for  this  rate  of  change  in  the  line  of 
sight  and  bring  it  to  zero,  the  missile  needs  to  turn  at  a 
corresponding  rate.  This  is  called  the  turn  rate  or  the 
acceleration  of  the  missile  and  is  given  by  the  following  equation: 

AM  =  SIGD  *  VM  *  n  (3.58) 

where 
n  =  Missile  navigation  constant 

E.    RELATIVE  COORDINATE  SYSTEM 
1.    Shift  of  Origin 

In  the  complex  and  randomly  occurring  scenario  of 
different  ship  headings,  missile  bearings  and  wind  directions, 
movement  of  the  chaff  after  full  bloom  can  take  many  paths  in  terms 
of  relative  motion.   In  order  to  minimize  the  disadvantageous 
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features  of  the  variables  involved,  we  have  decided  to  have  a 
relative  coordinate  system  centered  on  the  chaff  round.  This  is 
accomplished  with  the  help  of  the  following  equations: 

XR(h)  =  XS(h)  -  XC(h)  (3.59) 

YR(h)  =  YS(h)  -  YC(h)  (3.60) 

where  XR(N)  and  YR(N)  are  the  relative  abscissae  and  ordinates 

respectively.  A  sketch  depicting  the  shift  of  the  reference  system 

is  shown  in  Figure  16. 


^HIP 

(1.2) 
(-1.1) 
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CHAFF  (NEW  ORIGIN) 


(2.1) 
(0,0) 


(0,0) 
(-2,-1) 


Figure  16.  Shift  of  Reference  System 
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The  original  coordinates  of  the  positions  of  the  ship  and  chaff 
are  (1,2)  and  (2,1)  respectively  but  when  we  shift  the  origin  to 
the  chaff,  ship's  position  using  the  above  mentioned  equations 
becomes  (-1,1) . 

2.  Rotation  of  the  System 

The  next  step  is  to  rotate  the  new  coordinate  system  to 
the  bearing  of  the  missile  in  the  clockwise  direction.  Such  a 
rotation  is  desirable  and  handy  to  calculate  the  chaff  launch 
bearings  which  are  orthogonal  to  missile  direction.  This  is  done 
using  the  following  relationship: 

|     Cos  BRGM(h)         -Sin  BRGM(h)    | 

XRP(h)  =    |  |     XR(h) 

|     Sin  BRGM(h)  Cos  BRGM(h)     |    (3.61) 

where   XRP   and   XR   are   the   new  and  previous   state  vectors 

respectively.  Upon  multiplication  with  the  matrix  we  get  the 

following  equations: 

XRP(h)  =  XR(h)  *  Cos  BRGM(h)  -  YR(h)  *  Sin  BRGM(h)     (3.62) 

YRP(h)  =  XR(h)  *  Sin  BRGM(h)  +  YR(h)  *  Cos  BRGM(h)     (3.63) 

This  mathematical  model  can  be  better  understood  using  Figure  17. 

We  see  that  point  (0,2)  becomes  (-1,1.7)  when  the  system  is 

rotated. 

3 .  Relative  Velocity  Components 

The  next  important  step  is  to  calculate  the  relative 
velocity  components  using  the  initial  values  of  the  ship's  velocity 
and  the  velocity  of  the  wind  with  which  the  chaff  moves.  This  is 
shown  in  Figure  18. 
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MISSILE 


XRP 


Figure  17.  Rotated  Reference  System 
Using  Figure  18  we  get  the  following  eguations: 

XRPD(h)  =  VS*Sin(CSE(h)-BRGM(h) )  +  VW*Sin(GW  -  BRGM(h))  (3.64) 
YRPD(h)  =  VS*Cos(CSE(h)-BRGM(N) )  +  VW*Cos(GW  -  BRGM(h))  (3.65) 
where  VS  and  VW  are  the  velocities  of  ship  and  wind  respectively 
and  XRDP(h)/YRDP(h)  are  the  X  and  Y  components  of  the  relative 
velocity  respectively.  These  values  of  relative  velocity  components 
are  calculated  using  the  initial  conditions. 
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MISSILE 


WIND 


Figure  18.  Relative  Velocity  Components 

F.    PRE-COMPUTATION  METHOD 

The  random  variability  of  the  missile  direction,  ship's  course 
and  speed  at  the  time  of  missile  detection  along  with  wind  speed 
and  direction,  present  a  much  more  complex  problem  than  is  readily 
obvious.  It  is  extremely  difficult  to  obtain  the  best  solution  in 
all  the  possible  cases  with  conventional  chaff  alone.  Efforts  are, 
however,  made  in  this  work  to  enhance  optimization.  In  order  to 
achieve  this,  it  is  considered  useful  to  carry  out  a  rough  pre- 
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computation  which  can  give  an  idea  as  to  what  chaff  firing  bearing 
would  entail  the  largest  miss  distance.  Making  use  of  the  relative 
coordinate  system,  the  problem  becomes  somewhat  simpler  in 
formulation.  As  mentioned  earlier,  the  chaff  is  fired  orthogonal 
to  the  missile  bearing,  leaving  us  with  two  alternatives; one  to  the 
right  or  in  the  positive  XRP  direction  and  the  other  in  the 
negative  XRP  direction  i.e.,  to  the  left.  The  distance  traveled  by 
the  chaff  in  either  direction  is  computed  for  the  entire  run  of  the 
play.  A  couple  of  assumptions  made  in  this  regard  make  this  a 
rather  rough  but  a  pretty  good  estimate  of  the  situation.  The  first 
assumption  is  that  the  ship's  speed  is  taken  to  be  25  knots  which 
is  the  maximum  speed  of  the  ship  assumed  in  this  work  rather  than 
a  gradual  increase  from  the  ship's  initial  speed.  The  next 
assumption  ignores  the  alteration  to  the  first  leg  of  the  new 
courses  which  is  carried  out  to  escape  the  range  gate  in  minimal 
time.  A  direct  alteration  to  the  second  leg  from  the  ship's  initial 
course  is  assumed  instead,  which  is  perpendicular  to  the  missile 
direction.  The  following  equations  calculate  the  magnitude  of  the 
relative  motion  of  the  ship  with  respect  to  the  chaff  cloud  in 
either  direction. 

XRP1  =  XRP1  +  XRPD1*T  +  0.5*XRPDD*T2  (3.66) 

YRP1  =  YRP1  +  YRPD1*T  +  0.5*YRPDD*T2  (3.67) 

XRP2  =  -XRP1  +  XRPD1*T  +  0.5*XRPDD*T2  (3.68) 

YRP2  =  -YRP1  +  YRPD1*T  +  0.5*YRPDD*T2  (3.69) 

where 

XRPD1  =  XRPD1  +  XRPDD*T  (3.70) 


36 


YRPD1  =  YRPD1  +  YRPDD*T  (3.71) 

and 

XRPDD  =  -ACC*Sin(CSEDIFl)  (3.72) 

YRPDD  =  -ACC*Cos(CSEDIFl)  (3.73) 

where 

XRP1/YRP1  =  Distance  traveled  by  chaff  in  X/Y  direction 
XRPD1/YRPD1  =  Velocity  component  of  chaff  in  X/Y  direction 
XRPDD/YRPDD  =  Acceleration  component  of  chaff  in  X/Y  direction 
ACC  =  Acceleration  of  the  ship 

CSEDIF1  =  Difference  in  the  ship's  initial  course  and  the  new 
course 

This  is  shown  in  Figure  19. 

G.    RANGES  AND  BEARINGS 
1 .    Ranges 

There  are  three  different  ranges  of  interest  shown  in 
Figure  2  0  that  are  required  for  subsequent  calculations.  The  range 
between  the  missile  and  the  ship  center  is  given  by: 

R(h)  =  [  {XS(h)  -  XM(h)}2  +  {YS(h)  -  YM(h)}2  ]*  (3.73) 
Similarly  the  range  between  the  missile  and  the  chaff  center  is 
given  by: 

RCM(h)  =  [  (XC(h)  -  XM(h)  }2  +  (YC(h)  -  YM(h))2  ]*  (3.74) 
and  the  range  between  the  missile  and  the  centroid  or  the 
equivalent  target  is: 

RMT(h)  =  [  {XET(h)  -  XM(h)}2  +  (YET(h)  -  YM(h)}2  ]*     (3.75) 
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Figure  19.  Pre-Computed  Distances 
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Figure  20.  Various  Ranges  from  the  Missile 
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2 .    Bearings 

Like  the  ranges  mentioned  in  the  previous  section,  there 
are  also  some  bearings  that  are  necessary  for  subsequent 
calculations.  These  bearings  are  shown  in  Figure  21.  The  bearing 
of  the  equivalent  target  from  the  missile  is  given  by: 


BRGMT(h)  =  tan 


-1 


(h)  -  YM(h) 
(M  -  kM(h)  J 


YET (hi  -  YMi 
XEtT" 


(3.76) 
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Figure  21.  Various  Bearings  from  the  Missile 
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The  reciprocal  of  this  bearing,  i.e.,  the  bearing  of  the  missile 
from  the  equivalent  target  is  : 


BRGTM(h)  -  tan"   j^jnl  -   » 


The  bearing  of  the  ship  from  the  missile  is 


(3.77) 


BRGMS(h)    =   tan"1    < 


YS(h)    -   YM(h) 
Xsfhj    -   XM(hj 


(3.78) 


Lastly,  the  bearing  of  the  chaff  from  the  missile  is 


BRGMS(h)  =  tan' 


YSfh)  -  YM(h) 
XS(h)  -  XM(h) 


(3.79) 


H.    RANGE  AND  ANGLE  GATES 

The  range  gate  of  the  missile  is  approximately  a  function  of 
its  pulse  width  and  is  an  input  variable  to  the  problem.  The  angle 
gate,  on  the  other  hand,  is  dependent  upon  the  distance  of  the 
missile  from  the  target  and  the  half  power  beamwidth  of  the  missile 
radar  antenna.  This  is  calculated  using  the  following  equation: 

AG  =  RMT(h)  *  Tan(0B)  (3.80) 

We  also  need  to  calculate  the  distances  of  the  ship  and  the  chaff 
from  the  center  of  the  radar  resolution  cell,  which  is  the 
equivalent  target  or  the  centroid  of  the  two  radar  cross  sections. 
The  position  of  the  ship  with  respect  to  the  radar  resolution  cell 
center  along  the  range  and  angle  gate  directions,  is  calculated  as 
follows: 
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AGS  =  R(h)  *  Sin  |  BRGMS(h)  -  BRGMT(h)  |  (3.81) 

RGS  =  |  RMT(h)  -  R(h)  *  Cos{BRGMS(h)  -  BRGMT(h)}  |  (3.82) 

This  is  shown  in  Figure  22.  Similarly,  the  position  of  the  chaff 

is  given  by  the  following  equations: 

AGC  =  RMC(h)  *  Sin  |  BRGMC(h)  -  BRGMT(h)   |  (3.83) 

Resolution 
Cell 


Equivalent   Target 


Missile 
(XM,  YM) 

Figure  22.  Ship's  Position  Relative  to  Cell  Center 
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RGC  =  |  RMT(h)  -  RMC(h)  *  Cos{BRGMC(h)  -  BRGMT(h) }  | 
This  is  shown  in  Figure  23. 


(3.84) 


I.    RADAR  RESOLUTION  CELL  MODEL 

The  radar  resolution  cell  is  constructed  with  the  help  of  the 
range  and  angle  gates  of  the  missile  radar.  It  is  centered  on  the 
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Figure  23.  Chaff  Position  Relative  to  Cell  Center 
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equivalent  or  the  centroid  and  is  approximated  to  be  a  rectangle 
as  discussed  in  Chapter  II.  The  cell  is  oriented  towards  the 
missile  bearing.  As  mentioned  earlier,  the  range  gate  of  the 
missile  is  constant.  The  angle  gate,  on  the  other  hand,  is 
dependent  on  the  half  power  beamwidth  of  the  antenna  and  range  of 
the  missile  so  it  varies  as  the  missile  closes  in.  The  cell  is, 
therefore,  a  time  variant  rectangle  increasing  in  size  as  the 
missile's  range  decreases  and  reorienting  to  maintain  orthogonality 
with  the  missile  direction. 

The  depth  of  the  resolution  cell  is  twice  that  of  the  range 
gate  and  the  width  is  twice  that  of  the  angle  gate.  This  is  shown 
in  Figure  24. 
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Figure  24.  Radar  Resolution  Cell  Model  [Ref.  23] 
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The  edges  along  the  width  of  the  center  line  are  calculated 

as  follows: 

XW1  =  XET(h)  +  AG  *  Cos{ (BRGM(h) }  (3.85) 

YW1  =  YET(h)  -  AG  *  Sin{ (BRGM(h) }  (3.86) 

XW2  =  XET(h)  -  AG  *  Cos{ (BRGM(h) }  (3.87) 

YW2  =  YET(h)  +  AG  *  Sin{ (BRGM(h) }  (3.88) 

The  corners  of  the  cell  are  calculated  as  follows: 


XCELL(l)  =  XW1  +  RG  *  Sin{(BRGM(h) 

YCELL(l)  =  YW1  +  RG  *  Cos{(BRGM(h) 

XCELL(2)  =  XW1  -  RG  *  Sin{ (BRGM(h) 

YCELL(2)  =  YW1  -  RG  *  Cos{ (BRGM(h) 

XCELL(3)  =  XW2  -  RG  *  Sin{ (BRGM(h) 

YCELL(3)  =  YW2  -  RG  *  Cos{(BRGM(h) 

XCELL(4)  =  XW2  +  RG  *  Sin{(BRGM(h) 

YCELL(l)  =  YW2  +  RG  *  Cos{(BRGM(h) 


(3.89) 
(3.90) 
(3.91) 
(3.92) 
(3.93) 
(3.94) 
(3.95) 
(3.96) 


J.    PORTION  OF  THE  SHIP  INSIDE  THE  RESOLUTION  CELL 
1.    Ship  Fully  Inside  the  Cell 

The  entire  radar  cross  section  of  the  ship  lies  inside 
the  resolution  cell  initially.  After  firing  the  chaff,  the  ship 
attempts  to  maneuver  out  of  the  cell.  In  the  process,  the  center 
of  the  cell  shifts  to  the  equivalent  target.  The  new  situation 
alters  the  radar  cross  section  of  the  ship  as  seen  by  the  missile 
in  order  to  calculate  the  centroid.  This  cross  section  depends  on 
the  portion  of  the  ship  inside  the  cell  and  is  a  function  of  the 
following  parameters: 
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•  Radar  resolution  cell  center 

•  Size  of  the  cell 

•  Relative  location  of  the  ship  center 

•  Ship's  heading 

The  RCS  of  the  ship  is  assumed  to  be  the  maximum  presented  by 
the  ship  at  a  particular  aspect  angle  when  it  is  fully  within  the 
bounds  of  the  cell  as  shown  in  Figure  25. 
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Figure  25.  Limit  for  Maximum  RCS  of  the  Ship 
in  the  Resolution  Cell  [Ref.  23] 

2.    Ship  Partially  Inside  the  Cell 

When  some  part  of  the  ship  is  outside  the  cell,  we  have 
a  partial  radar  cross  section  case.  In  order  to  calculate  it,  we 
need  to  shift  our  reference  system  from  true  North  to  the 
horizontal.  This  shift  is  necessitated  by  the  fact  that  the 
computer  calculates  angles  from  the  horizontal  reference  as  shown 
in  Figure  26. 


45 


Inner  Circle 
With   Respect  to 
True  North 
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Outer  Circle 
With  Respect  to 
Fortran    Horizontal 
Reference 


+90. 


Figure  26.  Angle  Translation  to  Fortran  [Ref.  23] 

A  conversion  is,  therefore,  undertaken  to  cater  for  this  new 
reference  system.  We  define  a  new  angle  R,    as  shown  in  Figure  27 
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Figure  27.  True  North  to  Fortran  Conversion  [Ref.  23] 
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between  the  new  reference  and  the  direction  from  which  the  missile 
is  approaching  the  ship.  The  aspect  angle  of  the  ship  is  called 
'Saa'  and  it  is  calculated  using  8.  It  is  defined  as  the  angle 
between  the  ship's  heading  and  8  as  shown  in  Figure  28. 
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Figure  28.  Aspect  Angle  [Ref.  23] 

The  ship  is  approximated  as  a  rectangle  for  ease  of 
calculation.  Also  shown  in  this  figure  are  two  more  parameters; 
namely  Sir  and  S2r.  The  subscript  'r'  merely  represents  the  fact 
that  the  ship  is  exiting  through  the  range  gate  boundary  of  the 
cell.  If  it  exits  via  the  angle  gate  boundary,  then  we  replace  it 
with  subscript  'a'.  S2r  is  the  shortest  distance  between  the  range 
gate  boundary  and  the  ship  center  and  it  is  given  by: 

S2r  =  |  RG  -  RGS  |  (3.97) 

Sir  represents  the  portion  of  half  of  the  ship's  length  that  is 
inside  the  cell  and  is  calculated  with  the  help  of  the  following 
equation: 
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Slr  =  Cosflaa)  <3'98> 

In  the  case  shown  in  Figure  28,  where  the  ship  center  is 
inside  the  cell,  the  portion  of  the  ship's  RCS  in  the  cell  is  given 
by: 

PARTS  =  -J-  +  ffjfi  (3.99) 

where 

SI  =  Sir  or  Sla  depending  upon  how  the  ship  exits  the  cell 

Sarea  =  SL*SB 

In  case  the  ship  center  lies  outside  the  resolution  cell,  the 

calculation  is  as  follows: 


PARTS  =   4-  -  ffjfi  (3.100) 


If  the  ship  is  exiting  a  corner  (leaving  the  range  and  angle 
gate  boundaries  simultaneously) ,  then  the  amount  of  ship's  RCS  in 
the  cell  is  approximated  by  using  the  lesser  of  the  two  values 
obtained  from  Sir  or  Sla. 

K.    PORTION  OF  CHAFF  INSIDE  THE  CELL 
1.    Chaff  Fully  Inside  the  Cell 

The  chaff  is  fired  to  bloom  inside  the  radar  resolution 
cell.  After  it  blooms  the  wind  takes  over  and  the  chaff  cloud 
starts  moving  with  the  wind  vector.  Since  the  cell  is  centered  on 


48 


the  equivalent  target,  some  portion  of  the  chaff  might,  dependent 
on  the  situation,  spill  out  of  the  cell.  In  order  to  get  an 
estimate  of  the  location  of  the  center  of  the  cell,  it  is  important 
to  know  how  much  of  chaff  is  inside  the  cell.  This  depends  on  the 
following  parameters: 

•  Relative  location  of  chaff  center 

•  Size  of  chaff 

•  Location  of  resolution  cell  center 

•  Size  of  the  cell 

The  chaff  lies  completely  within  the  cell  when  the  distance 
of  its  center  from  the  cell  bounds  is  at  least  equal  to  the  radius 
of  the  chaff.  This  is  shown  in  Figure  29. 


CRAD 


Chaff  Cloud 


Figure  29.  Limit  of  Maximum  Chaff  RCS  in  the  Cell  [Ref.  23] 
2.    Chaff  Partially  Inside  the  Cell: 

We   define   the   following   parameters   to   assist   in 
calculation  of  the  chaff's  RCS  inside  the  cell: 


Carea  =  r    *    (CRAD) 


(3.101) 
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where 

CRAD  =  Chaff  radius 

Rlr  =  |  RGC  -  RG  |  (3.102) 

Rla  =  |  AGC  -  AG  |  (3.103) 

If  the  chaff  is  located  at  the  range  gate  boundary,  the 

situation  will  be  as  shown  in  Figure  30.  In  this  particular  case, 

more  than  half  the  chaff  lies  outside  the  bounds  of  the  cell. 


Chaff 


(XC,  YC) 


Figure  30.  Portion  of  Chaff  RCS  in  the  Cell  [Ref.  23] 

If  we  expand  the  view  of  interest,  as  shown  in  Figure  31,  we  can 
calculate  a  few  more  parameters  that  would  assist  us  in  finding 
chaff  RCS  inside  the  cell. 
We  see  that: 


A2r  =  Cos 


-1 


Rlr 

CRAD 


(3.104) 
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The  sector  bounded  by  twice  this  angle  is  given  by: 
SECTr  =  1/2  (2*A2r)  CRAD2 


(3.105) 
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Resolution   Cell 
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Figure  31.  Partial  Chaff  RCS  in  the  Cell  Blown  Up  [Ref.  23] 

Now 

R2r  =  CRAD  *  Sin(A2r)  (3.106) 

The  area  of  the  two  triangles,  which  is  also  the  portion  of  the 
sector  outside  the  radar  resolution  cell  is: 

TRIr  =  2(I/2*Rlr*R2r)  =  Rlr  *  R2r  (3.107) 

In  case  the  chaff  exits  the  cell  through  the  angle  gate, 
the  same  procedure  is  applied  using  subscript  Rla.  If  both  the 
gates  are  cut  off  simultaneously,  then  we  use  both  procedures  to 
calculate  the  area  or  the  chaff  RCS  inside  the  cell.  All  three 
possible  cases  are  shown  in  Figure  32. 

When  Rlr  or  Rla  are  zero  individually,  the  portion  inside 
the  cell  is  one  half  of  the  total.  However,  when  they  are  both  zero 
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simultaneously,  as  in  the  corner  case,  then  it  is  only  one  quarter 
of  the  total. 

PARTc  =  0.5 


Chaff 


PARTc  =  0.25 


Figure  32.  Various  Cases  of  Partial  Chaff  RCS  [Ref.  23] 

If  RG  >RGC,  then  most  of  the  chaff  cloud  lies  inside  the 
resolution  cell.  The  portion  of  chaff  in  the  cell  is  calculated  as 
follows: 


PARTC  = 


CAREA  -  SECTr  +  TRIr 
CAREA 


(3.108) 


IF  RG  <  RGC,  then  the  RCS  in  the  cell  is  less  than  half  and  is 
given  by: 


partp  -   SECTr  -TRIr 
FAK1C        CAREA 


(3.109) 


When  the  chaff  exits  the  cell  from  a  corner  we  have  four 
situations: 

1.  RGC  <  RG  and  AGC  <AG 
This  is  shown  in  the  Figure  33  along  with  the  actual  and 
approximated  radar  cross  sections. 


52 


Angle 
Gate 


Angle 
Gate 


RCS 
Approximation 


Actual  RCS  within  the  Cell 


Actual  RCS 
within 
the  Cell 


Range 
Gate 


Range 
Gate 


Approximation 


Figure   33.    Chaff   Exiting  the   Cell   Corner    [Ref.    23] 

2 .    RGC    <    RG   and   AGC    >   AG 
This   situation    is    shown    in   Figure    34. 
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Figure    34.    RGC.LT.RG    &   AGC.GT.AG    [Ref.    23] 
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3.  RGC>  RG  and  AGC  <  AG.  This  is  shown  in  Figure  35. 
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Figure  35.  RGC.GT.RG  &  AGC.LT.AG  [Ref.  23] 

4 .  RGC  >  RG  and  AGC  >  AG 
In  this  case  the  larger  of  the  two  areas  produced  in  the  cell  by 
RGC  or  AGC  is  approximated  to  be  the  chaff  RCS  of  interest  as  shown 
in  Figure  36.  Less  than  a  guarter  of  the  chaff  is  inside  the  cell 
and  it  is  taken  to  be  one  half  the  area  cut  off  by  the  lesser  AGC 
or  RGC. 
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Figure    36.    RGC.GT.RG   &   AGC.GT.AG    [Ref.    23] 
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L.    CHAFF  LAUNCH  BEARING  AND  SHIP  MANEUVERS 

1.  General  Description 

In  order  to  maximize  time  to  break  lock  and  miss 
distance,  a  double  maneuver  by  the  target  ship  is  simulated.  The 
first  maneuver  is  to  get  out  of  the  radar  resolution  cell  of  the 
missile  via  the  shortest  route,  i.e,  either  towards  the  missile  or 
directly  away  from  it.  The  next  maneuver  is  to  alter  to  a  course 
perpendicular  to  the  missile  direction  when  break  lock  is  achieved 
and  the  ship  is  outside  the  resolution  cell  of  the  missile's  radar. 

2 .  Chaff  Launch  Bearing 

First  of  all,  the  bearing  to  launch  the  chaff  is  worked 
out.  It  is  basically  desired  to  launch  the  chaff  at  right  angles 
to  the  missile  direction.  The  two  possible  true  bearings  are 
calculated  as  follows: 

BRGC1  =  BRGM(l)  +  90  (3.110) 

BRGC2  -  BRGM(l)  -  90  (3.111) 

The  decision  regarding  the  chaff  firing  bearing  is 
dependent  upon  the  results  of  the  pre-computation  in  order  to 
optimize  the  miss  distance.  If  the  ship  travels  more  along  the 
positive  XRP  direction,  we  fire  the  chaff  towards  BRGC1  and  vice 
versa. 

In  the  particular  case  shown  in  Figure  37,  BRGC  is  equal 
to  BRGC2 .  It  must  be  kept  in  mind  that  the  missile  is  shown  to  come 
from  the  top  because  of  the  rotated  reference  system  and  could  in 
actuality  be  approaching  from  any  direction  in  the  azimuth. 
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Figure  37.  Pre-Computation  for  Calculating  BRGC 
3.    First  Maneuver 

The  next  object  is  to  calculate  the  first  course  to  steer 
to  get  out  of  the  radar  resolution  cell  after  having  fired  the 
chaff.  This  situation  is  shown  in  Figure  38. 

The  choice  to  select  one  of  the  course  options  is  based 
on  the  relative  movement  of  the  ship  along  the  vertical  axis,  i.e., 
YRP  direction  which  is  oriented  along  the  missile  line  of  flight. 
The  following  algorithm  in  the  program  enables  us  to  decide  what 
course  option  to  take: 
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missile. 


If  YRP(N)  >  0  then  CSE2  or  the  first  leg  is  into  the 


If  YRP(N)  <  0  then  CSE2  is  away  from  the  missile. 
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Figure  38.  Course  to  Steer 
4 .    Second  Maneuver 

Once  the  ship  is  out  of  the  radar  resolution  cell  we 
calculate  the  next  course  to  steer  to  maximize  miss  distance.  This 
is  initiated  when  one  of  the  following  conditions  occur: 

a)  YRP(N)  >  RG 

b)  XRP(N)  >  AG 

In  either  case,  the  object  is  to  ascertain  whether  the 
ship  is  out  of  the  radar  resolution  cell,  implying  that  break  lock 
has  occurred.  However,  the  subroutine  for  calculating  the 
parameters  for  the  second  maneuver  only  calculates  CSE3  and  the 
turn  direction.  The  actual  implementation  of  affecting  this  course 
alteration  based  on  the  two  conditions  already  mentioned,  is 
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achieved  in  the  main  program.  Since  the  second  leg  of  the  new 
course,  CSE3,  is  at  right  angles  to  the  first  leg,  CSE2,  there  are 
again  two  choices  presented  as  shown  in  Figure  39. 

It  has  been  supposed  here  that  CSE2  is  equal  to  CSE0PT2, 
referring  to  Figure  38. 
Moreover, 

CSE0PT(3)  =  BRGM(l)  +  90  (3.112) 

CSE0PT(4)  =  BRGM(l)  -  90  (3.113) 
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Figure  39.  Second  Leg  of  the  Course 

The  selection  of  CSE3  from  these  two  options  is  related 
to  the  bearing  in  which  the  chaff  is  launched  and  is  usually  in  its 
opposite  direction.  In  the  situation  shown  in  Figure  39, 

CSE3  =  CSE0PT(4) 

The  direction  to  turn  is  calculated  in  a  manner  similar 
to  the  one  described  in  Section  3. 
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IV.   DESCRIPTION  OF  PROGRAM 

The  simulation  program  is  written  in  FORTRAN-77  (Microsoft 
version  3.2).  In  order  to  show  the  run  of  play  at  various  times, 
the  program  uses  a  compatible  Plot-88  (1984)  graphics  package. 
A.    INITIAL  CONDITIONS 

The  initial  conditions  are  set  up  through  an  interactive 
exchange  with  the  operator  in  order  to  input  all  the  important 
variables.  The  following  variables  are  used  as  inputs  and  they  have 
been  assigned  values  that  are  considered  reasonably  standard  but 
they  can  be  altered  through  a  simple  modification  in  the  program. 

TABLE  4 
INPUT  VARIABLES  WITH  TYPICAL  VALUES 


S.No 

Variable 

Description 

Value 

1. 

CRAD 

Radius  of  Chaff 

50  yards 

2. 

SL 

Ship  Length 

200  mtrs 

3. 

SB 

Ship  Breadth 

20  meter 

4. 

VSMAX 

Max  Ship  Velocity 

2  5  knots 

5. 

N 

Missile  Nav  Constant 

3 

6. 

VS(1) 

Initial  Ship  Vel 

2  0  knots 

7. 

VMM 

Missile  Velocity 

600  knot 

8. 

TR 

Turn  Rate  of  the  Ship 

5'/s 

9. 

VR 

Ship's  Acceleration 

0.5  kts/s 

10. 

AM 

Missile  Acceleration 

30  yards 

11. 

CDIST 

Range  of  Chaff 

200  yards 
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ITTG  =  Integer  i   TTG  +  f10^}00) 


(4.2) 


C.    MAIN  BODY  OF  THE  PROGRAM 

The  main  body  of  the  program  iterates  through  an  array  with 
a  maximum  of  60  iterations  allowed.  The  values  of  missile,  chaff 
and  ship  parameters  are  updated  along  with  other  related  parameters 
in  each  time  step.  The  related  parameters  include  gates,  angles, 
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A  scaling  factor  is  included  in  the  input  variables  to  enable 
actual  ranges,  velocities  and  accelerations  to  be  depicted  on  the 
screen  of  the  monitor.  The  screen  center  (4.0,3.0)  is  used  as  a 
reference  with  (0.0,0.0)  being  the  lower  left  hand  corner  and 
(8.0,6.0)  the  upper  right  hand  corner  of  the  screen. 

B.    RUN  TIME 

The  program  is  incremented  in  steps  of  one  second  intervals 
until  the  missile  gets  to  within  500  yards  from  the  ship. 
Thereafter,  to  show  details  of  action,  we  use  0.1  second  intervals. 
The  total  time  to  go  is  found  by  the  following  relationship: 

TTG  =   B^E  (4#1) 

where 

RANGE  =  Initial  Range  of  the  Missile 

The  total  program  step  intervals  including  0.1  second  steps  are 

given  by: 


equivalent  target  parameters  and  radar  cross  sections,  etc.  The 
missile  radar  seeker  is  turned  off  when  its  range  from  the 
equivalent  target  is  250  yards  resulting  in  zero  missile 
acceleration  and  from  thereon,  the  missile  motion  is  essentially 
ballistic.  When  the  missile  has  flown  through  the  chaff  cloud,  its 
range  starts  increasing  again  and  the  main  loop  is  completed  on 
that  iteration.  If  the  missile  passes  within  50  yards  of  the  ship, 
it  is  considered  to  have  hit  the  ship.  The  implementation  of 
alteration  of  courses  to  maneuver  the  ship  is  also  achieved  in  the 
main  body  of  the  program  using  the  decision  algorithms  mentioned 
in  Section  '  L' . 

D.    USER  INTERACTION  AND  OUTPUT 

After  having  initialized  all  the  variables,  the  user  is 
presented  with  two  options;  either  to  let  the  program  calculate  the 
control  parameters  by  resorting  to  auto  mode  or  to  manually  inject 
these  parameters.  The  parameters  consist  of  chaff  firing  bearing, 
both  legs  of  the  courses  to  be  steered  to  maneuver  the  ship  and  the 
turn  directions.  Subsequent  to  this  stage,  the  user  is  presented 
with  a  picture  comprising  of  initial  conditions  as  shown  in  Figure 
40.  The  essential  parameters  are  listed  on  the  picture  to  its 
right  and  at  the  bottom.  Thereafter,  the  user  is  presented  with 
another  menu  of  the  following  options: 

1.    Option  One 

This  option  draws  the  missile,  ship  and  chaff  paths  and 
gives  the  miss  distance  at  the  bottom  as  shown  in  Figure  41. 
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Figure   40.    Initial   Conditions 
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Figure  41.  Paths  Traversed  by  Missile,  Ship  &  Chaff 
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2.    Option  Two 

This  option  provides  a  pictorial  representation  quite 
similar  to  the  initial  conditions  picture  and  we  see  the  situation 
in  the  resolution  cell  in  5  second  increments.  It  depicts  the  run 
of  play  showing  successive  positions  of  the  missile,  chaff  and  the 
ship.  An  example  is  shown  in  Figure  42. 
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Figure  42.  Instantaneous  Position  of  the 
Resolution  Cell 
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3.    Option  Three 

Option  three  facilitates  changes  in  the  input  scenario 
parameters  with  the  help  of  a  menu.  It  is  incorporated  to  provide 
flexibility  to  the  user  to  make  any  alterations  in  the  input 
parameters  as  desired.  The  parameters  can  be  changed  by  choice 
either  individually  or  in  totality.  The  menu  of  options  is  as 
follows: 

TABLE  5 
CHANGEABLE  PARAMETERS 


MISSILE  BEARING 

MISSILE  RANGE 

SHIP'S  HEAD 

SCREEN  SCALE 

MISSILE  RADAR  BEAMWIDTH 

RANGE  GATE 

MAXIMUM  SHIP  RADAR  CROSS  SECTION 

MINIMUM  SHIP  RADAR  CROSS  SECTION 

RADAR  CROSS  SECTION  OF  CHAFF 

TRUE  WIND  DIRECTION 

TRUE  WIND  SPEED 

ALL  PARAMETERS 

RETURN  TO  MAIN  PROGRAM 


4 .    Option  Four 

This  option  terminates  the  program. 

E.    SUBROUTINES 

A  number  of  subroutines  are  made  use  of  to  break  the  program 
into  smaller  modules.  Their  respective  functions  and  broad  outlines 
are  presented  here  to  give  an  overview  to  the  reader. 
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1.  Auto  Subroutines 

There  are  two  auto  subroutines.  The  function  of  the  first 
auto  subroutine  is  to  calculate  the  chaff  launch  bearing,  BRGC,  the 
first  leg  of  the  ship  maneuvers,  CSE2  and  the  direction  to  enable 
the  shorter  turn  from  the  initial  course.  In  order  to  optimize  miss 
distance,  a  fairly  accurate  estimate  of  the  chaff  launch  bearing 
is  obtained  using  the  pre-computation  method.  The  decision  related 
to  CSE2  calculation  is  based  on  the  algorithm  in  section  'L'.  The 
second  subroutine  works  out  the  second  leg  of  the  maneuver  and  the 
direction  to  turn  away  from  the  chaff. 

2.  Ship,  Rescell  And  Chaff  Subroutines 

These  subroutines  calculate  the  coordinates  of  the 
corners  of  the  resolution  cell,  the  hexagon  representing  the  chaff 
cloud  and  the  polygon  representing  the  ship. 

3.  Shipcell  And  Chaff cell  Subroutines 

These  subroutines  calculate  the  radar  cross  sections  of 
the  ship  and  the  chaff  inside  the  resolution  cell  to  enable 
calculation  of  the  centroid  in  the  main  program. 

4 .  Aspect  Subroutine 

The  aspect  angle  of  the  ship  as  it  exits  the  resolution 
cell  is  passed  on  to  the  Shipcell  subroutine  for  subsequent  RCS 
calculations. 

5.  Missile,  Wind,  Path  And  Table  1  Subroutines 

These  are  PLOT-88  subroutines.  The  first  three  are  used 
to  draw  the  missile  line,  wind  line  and  paths  traced  by  the 
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missile,  ship  and  chaff.  Table  1  subroutine  constructs  a  table  of 
essential  parameters  when  Screnl  &  2  subroutines  are  used. 
6.   Scren  1  2  &  3  Subroutines 

Scren  1  subroutine  draws  the  missile  line  of  flight,  wind 
vector  and  the  ship  polygon  based  on  the  initial  conditions  input 
by  the  operator.  Scren  2  subroutine  is  used  to  draw  the  same 
picture  but  at  5  second  intervals  to  offer  a  look  at  the  situation 
dynamically.  Scren  3  subroutine  enables  the  operator  to  change 
desired  parameters  and  updates  the  memory  with  their  latest 
version. 
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V.   RESULTS/RECOMMENDATIONS 

A.    RESULTS 

The  random  variability  of  the  parameters  associated  with  the 
ship,  missile  and  the  chaff  in  a  given  battle  scenario  can  result 
in  a  very  large  number  of  possible  situations  at  sea.  Some 
combinations  of  the  parameters  mentioned  can  produce  a  more 
favorable  situation  for  the  ship  than  others  to  thwart  the  anti- 
ship  missile.  In  some  cases,  a  single  chaff  round  deployed  by 
itself  would  have  reduced  effectiveness  or  might  not  be  effective 
at  all.  A  strong  wind  blowing  from  the  direction  of  the  missile  or 
from  the  opposite  side  will  move  the  chaff  via  the  shortest  path 
out  of  the  resolution  cell.  In  such  cases,  the  chaff  will  exit  the 
cell  before  the  ship  can  do  so;  thus  foiling  the  effort  to  break 
lock.  Firing  the  chaff  into  the  predicted  range  gate  would  avoid 
this  problem.  It  is  not  considered  appropriate  to  discuss  results 
produced  by  the  algorithm  used  in  this  work  for  all  the  possible 
cases.  Some  random  situations  have,  therefore,  been  selected  to 
give  the  reader  an  idea  of  the  working  of  the  model  developed. 

1.    Situation  1 

The  salient  input  variables  for  this  situation  are  as 
follows: 

•  Missile  Bearing  -  000  Degrees 

•  Missile  Range    -  10,000  Yards 

•  Ship's  Course    -  225  Degrees 

•  Wind  Direction   -  330  Degrees 

•  Wind  Speed      -  15  Knots 
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The  initial  picture  is  shown  in  Figure  43.  The  bearing  to  launch 
as  calculated  by  the  program  when  selected  in  automatic  mode  is 
090.  This  has  been  worked  out  by  the  pre-computation  technique.  If 
the  manual  mode  is  selected  instead,  the  user  is  presented  the  same 
picture  without  the  launch  bearing  and  asked  to  input  the 
following: 

•  Bearing  to  launch  the  chaff 

•  First  course  to  steer  to  exit  the  range  gate 

•  Direction  to  turn  to  first  course 

•  Second  course  to  steer  to  get  away  from  the  missile  line  of 
fire 

•  Direction  to  turn  to  second  course 

In  the  automatic  mode,  however,  the  program  works  out  the 
solution  itself  and  then  the  result  of  the  simulation  can  be  viewed 
in  two  ways.  The  first  graphic  display  shows  the  paths  traced  by 
the  ship,  the  missile  and  the  chaff.  This  is  shown  in  Figure  44. 
In  the  case  shown,  the  missile  has  missed  the  ship  by  882.9  yards. 
The  ship  first  maneuvered  to  steer  southwards  to  break  lock  ahd 
then  westwards  to  maximize  miss  distance.  The  second  mode  of 
presentation  displays  the  situation  every  5  seconds  until  the 
missile  hits  one  of  the  two  possible  targets.  The  situation  after 
15  seconds  is  shown  in  Figure  45. 
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Figure  43.  Initial  Conditions (Situation  1) 
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Figure  44.  Paths  display (Situation  1) 
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Figure  45.  Instantaneous  Positions (Situation  1) 
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2.    Situation  2 

The  input  variables  are  as  follows: 

•  Missile  Bearing  -  280  Degrees 

•  Missile  Range   -  10,000  Yards 

•  Ship's  Course   -  150  Degrees 

•  Wind  Direction   -  220  Degrees 

•  Wind  Speed      -  15  Knots 

Figure  46  shows  the  initial  conditions.  The  bearing  to  launch  the 
chaff  is  010.  The  paths  traversed  and  the  situation  after  15 
seconds  are  shown  in  Figures  47  and  48  respectively. 
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Figure  46.  Initial  Conditions (Situation  2) 
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Figure  47.  Paths  display (Situation  2) 
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Figure  48.  Instantaneous  Positions (Situation  2) 
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B.    RECOMMENDATIONS 

The  deployment  of  chaff  in  seduction  mode  poses  a  complicated 
problem  in  terms  of  the  randomness  of  those  variables  that  are  not 
controllable  by  the  ship.  They  include  the  missile's  bearing,  the 
beam  width  of  its  radar,  the  ship's  course  at  the  time  of  initial 
detection  and  the  wind  vector.  Combined  with  the  initial  detection 
range  of  the  missile  these  factors  determine  the  size  of  the  radar 
resolution  cell  and  as  a  result,  the  time  spent  by  the  chaff  round 
inside  the  cell.  This  is  also  the  time  duration  available  to  the 
ship  to  vacate  the  cell  and  is,  therefore,  critical  in  solving  the 
problem. 

The  variables  that  can  be  controlled  by  the  ship  are  the  chaff 
range  and  bearing,  ship's  course  and  speed,  its  turn  rates  and 
acceleration.  The  last  three  have  a  limited  contribution  in  this 
context.  The  control  of  the  ship's  heading  and  the  bearing  to  fire 
the  chaff  is  of  the  utmost  significance.  In  order  to  increase  the 
effectiveness  of  the  control  variables  and  to  make  the  algorithm 
developed  more  adaptable,  some  recommendations  are  made  in  the 
ensuing  paragraphs. 

1.    Multiple  Chaff  Firing 

In  order  to  tackle  the  cases  where  the  chaff  exits  the 
radar  resolution  cell  before  the  ship,  multiple  chaff  rounds  must 
be  launched.  The  second  chaff  should  be  reseeded  inside  the  cell 
prior  to  the  exit  of  the  first.  Its  initial  placement  should  depend 
on  how  far  the  ship  is  from  the  cell  bounds  when  the  first  round 
is  about  to  exit  or  how  much  time  it  requires  at  the  point  to  leave 
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the  cell  keeping  the  bloom  time  delay  in  mind.   Thereafter, 
subsequent  chaff  rounds  may  be  fired  in  the  same  way  as  required 
until  the  ship  breaks  lock. 
2.    Chaff  Walk  Off 

This  is  another  philosophy  involving  multiple  chaff 
launches.  If  multiple  chaff  rounds  are  fired  on  the  same  bearing 
but  at  increasing  ranges,  rapid  break  lock  could  be  achieved  and 
then  the  ship  would  need  to  carry  out  only  one  maneuver  to  increase 
miss  distance.  The  chaff  rounds  can  be  fired  simultaneously  or  in 
quick  succession  so  that  they  are  waiting  in  anticipation  and  the 
centroid  shifts  in  a  discrete  manner.  To  explain  this,  let  us  study 
an  example  shown  in  Figure  49. 

With  the  missile  approaching  from  the  north,  the  first 
chaff  is  placed  at  the  cell's  angle  gate  boundary  on  a  bearing 
orthogonal  to  missile  direction.  The  bearing  can  be  offset  to  cater 
for  the  bloom  time  delay.  If  the  chaff  radar  cross  section  is  equal 
to  that  of  the  ship,  the  centroid  will  be  established  at 
approximately  the  mid-point  of  the  two,  called  XI.  The  resolution 
cell  now  moves  to  the  right  and  is  called  Cell  'B'.  At  this  point 
it  includes  the  ship,  Chaff  1  and  Chaff  2.  The  centroid  now  shifts 
to  X2  and  the  cell  to  Cell  'C.  The  third  chaff  round  displaces 
the  cell  position  to  Cell  'D'  and  the  centroid  to  X3 .  It  can  be 
seen  that  the  ship  is  now  very  close  to  the  left  hand  angle  gate 
boundary  and  can  steer  out  of  the  cell  to  achieve  a  rapid  break  and 
then  increase  the  miss  distance.  We  feel  that  this  technique  can 
enhance  the  effectiveness  of  chaff  in  anti-ship  missile  defense. 
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Figure  49.  Chaff  Walk  Off 

3.  Multiple  Solutions 

In  order  to  achieve  optimum  results,  the  program  can  be 
used  in  manual  mode  to  inject  different  control  inputs.  The  program 
run  time  is  quite  short  and  a  slight  modification  to  compare 
results  of  both  automatic  and  manual  modes  can  recommend  the  best 
solution  to  be  adopted  in  a  particular  situation. 

4.  Battle  Group  Operation 

When  operating  in  a  battle  group,  non-judicious 
deployment  of  chaff  can  result  in  seducing  the  chaff  onto  a 
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friendly  unit.  It  is,  therefore,  imperative  to  introduce  safe 
firing  bearings  in  the  algorithm  to  ensure  safety  of  ships  in 
company . 

5.  Fixed  Chaff  Launcher 

The  program  needs  to  be  adapted  to  include  units  equipped 
with  fixed  chaff  launchers  to  optimize  firing  bearings  while 
minimizing  maneuvers  by  the  ship.  The  complexity  of  the  problem  is, 
of  course,  increased  and  the  ship's  actions  are  further  limited. 
If  there  are  limitations  of  CIWS  firing  bearings,  they  can  also  be 
included  to  allow  for  an  integrated  response. 

6.  Integration  with  Active  ECM 

If  the  ship  is  equipped  with  an  active  ECM  capability, 
integration  of  chaff  deployment  with  the  active  ECM  response  to  an 
anti-ship  missile  threat  can  be  a  decisive  factor.  The  break  lock 
may  first  be  achieved  with  active  ECM  and  chaff  can  then  be  placed 
between  the  ship  and  the  false  target  to  deceive  the  missile  and 
also  to  increase  miss  distance.  Active  decoys  such  as  Nulka  and 
Winnin  employ  similar  techniques. 

7.  Stochastic  Modelling 

Due  to  the  randomness  of  the  variables  involved  in  the 
anti-ship  defense  problem,  stochastic  modelling  may  be  employed  to 
achieve  optimum  results.  In  this  paper,  wind  direction  and  speed 
are  assumed  to  be  constant.  This  is  usually  not  true  in  real  sea 
conditions.  Noise  in  the  wind  vector  can  also  be  introduced  through 
an  appropriate  probability  distribution  to  depict  more  realistic 
environmental  conditions. 
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8.  Radar  Cross  Section  Calculations 

The  mathematical  modelling  to  calculate  radar  cross 
sections  of  both  the  ship  and  the  chaff,  use  some  approximations 
when  a  simultaneous  exit  from  the  angle  and  range  gates  takes 
place.  Also,  in  case  of  simultaneous  spill-over  from  opposing  sides 
of  the  cell  through  the  range  gate,  the  partial  radar  cross  section 
is  calculated  assuming  that  the  spill-over  is  from  one  side  only. 
The  cross  section  which  lies  outside  the  cell  on  the  opposite  side 
is  thus  neglected.  The  reader  should  be  aware  of  these 
approximations  and  make  suitable  changes  if  so  desired. 

9.  Programming  Modifications 

The  simulation  program  is  written  in  FORTRAN  which  has 
some  inherent  limitations.  It  is  recommended  that  a  more  powerful 
and  flexible  language  be  used  to  improve  upon  the  program 
efficiency  and  to  incorporate  the  changes  suggested  more 
effectively.  The  dimension  figure  of  the  arrays  can  be  increased 
to  include  larger  initial  detection  ranges  which  may  be  applicable. 
Parameters  such  as  chaff  radius,  ship's  length,  turning  rates  and 
acceleration/deceleration  tables  may  be  included  as  input  variables 
to  adapt  the  program  to  given  chaff  and  ship  charateristics. 

C.    CONCLUSIONS 

Resurgence  of  decoys  in  recent  years  has  highlighted  their 
importance  in  the  area  of  anti-ship  missile  defense.  Although  the 
missile  radars  are  becoming  smarter  by  including  advanced  ECCM 
features,  chaff  is  considered  to  be  a  viable  defensive  technique 
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even  if  deployed  in  a  stand  alone  manner  such  as  centroid  mode. 
This  is  especially  valid  for  the  smaller  navies  who  have  limited 
resources.  Improvement  in  chaff  deployment  algorithms  and,  when 
possible,  integrated  operations  with  active  ECM  and  CIWS  assets  can 
make  chaff  a  potent  defensive  measure.  It  is  simple  to  use  and 
relatively  inexpensive  among  the  known  type  of  decoys.  Its 
intelligent  use  can  make  it  effective  and  dependable.  Although  wars 
have  not  been  won  by  a  single  factor  alone,  chaff  can  certainly 
become  one  of  the  important  factors  in  a  limited  conventional 
conflict. 
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APPENDIX.  SIMULATION  PROGRAM 
C 

C  ******************************************* 

C  *      DEPLOYMENT  OF  CHAFF  IN  CENTROID      * 

C  *      MODE  AGAINST  ANTI  SHIP  MISSILES      * 

C  *  USING  A  * 

C  *       VARIABLE  AZIMUTH  AND  ELEVATION      * 

C  *  LAUNCHER  * 

C  *  24  AUG  89     SIDDIQ.FOR  * 

C  *  M.  SIDDIQ  * 

C  ******************************************* 

c 
c 
c 
c 

C  THE  PROGRAM  HAS  NO  PRINTER/  PLOTTER  OPTIONS 

C  ITS  OUTPUT  GOES  TO  THE  MONITOR  ONLY 

C  IT  HAS  AN  AUTOMATIC  OR  MANUAL  CHAFF  LAUNCHING  OPTION 

C  IT  CREATES   NO  DATA  FILE 

C 

C    VARIABLES  DECLARATION 

c    


COMMON/ TIME/ I TTG , TTG 
COMMON/ TME/T 

COMMON/TRACK/  XM(60) ,YM(60) ,XS(60) ,YS(60) ,XC(60) ,YC(60) 
COMMON/WD/ GW , VW 
COMMON/MISS/RG , THETA , AG , PI 

COMMON/ABCD/BRGM(60) , R( 60) , SC, VS (60) ,CSE(60) ,BRGMT(60) , 
*BRGTM(60) ,RMT(60) 
COMMON/EQ/  XET(60) ,YET(60) 

COMMON/ SHCL/ AGS , RGS , BRGMS ( 60 ) , SAA , PARTS , SAAN ( 60 ) 
COMMON/ CFCL/ PARTC , CRAD , RGC , AGC 
COMMON/ SHP/ S  L , S  B 

COMMON/ RC/RCSS , RCSSP ( 60 ) , RCSC , RCSCP (60) 
COMMON/MINN/RMIN 
COMMON/ CHF/BRGC, I SHOOT 
COMMON/ TURN/ TR, CSEOFF 
COMMON/ FG/ I FLAG 1 (60) , IFLAG2 (60) 
COMMON/ PRN/ I OPORT, MODEL 

COMMON/DOTS/XMD(60) ,YMD(60) ,XMDD(60) ,YMDD(60) ,XSD(60) ,YSD(60) , 
*XSDD(60) ,YSDD(60) 

COMMON/REL/XR(60) ,YR(60) ,XRD(60) ,YRD(60) ,XRP(60) ,YRP(60) ,XRPD(60) , 
*YRPD(60) 
REAL  AMM,NAV,T 
REAL  CDIST 
REAL  TTG,GT,X,Y 
REAL  SIG,SIGD,AM,GM 
IOPORT  =  99 
MODEL  =99 
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c 

c 

c 

MISSILE: 

c 

c 

BRGM 

c 

BRGMS 

c 

BRGMT 

MISSILE (REF) 

C 

BRGMC 

C 

R/R(N) 

C 

RANGE 

C 

R(1)/RMIN 

C 

RMT(N) 

C 

RMC 

C 

VMM 

C 

VM 

C 

GM 

C 

AM 

C 

AMM 

C 

XM 

C 

YM 

C 

XMD 

c 

YMD 

c 

XMDD 

c 

YMD 

c 

SHIP: 

c 

c 

SL 

c 

SB 

c 

CSE 

c 

CSE2 

c 

CSE3 

c 

CSEOPT 

c 

CSEDIF1 

c 

TOPT 

c 

GS 

c 

VS 

c 

AS(VR) 

c 

TR 

c 

U(ALSO  Ul 

c 

SAA 

c 

XS 

c 

YS 

c 

XSD 

c 

YSD 

c 

XSDD 

c 

YSSD 

c 

RCSS 

c 

RCSSMAX 

c 

RCSSMIN 

c 

RCSSP(l) 

c 

RCSSP(N) 

LIST  OF  VARIABLES 


BEARING  OF  MISSILE  FROM  SHIP(WRT  NORTH) 
BEARING  OF  SHIP  FROM  MISSILE (WRT  REFERENCE) 
BEARING  OF  EQUIVALENT  TARGET  FROM 

BEARING  OF  CHAFF  FROM  MISSILE (REF) 

INSTANTANEOUS  MISSILE  RANGE  FROM  SHIP 

INITIAL  MISSILE  RANGE 

SCALED  INITIAL  MISSILE  RANGE 

MISSILE  RANGE  FROM  EQUIVALENT  TARGET 

MISSILE  RANGE  FROM  CHAFF 

MISSILE  VELOCITY 

SCALED  MISSILE  VELOCITY 

MISSILE  DIRECTION 

MISSILE  ACCELERATION 

SCALED  MISSILE  ACCELERATION ( 3 O/SC) 

X-COORD  OF  MISSILE  POSITION 

Y-COORD  OF  MISSILE  POSITION 

X-COORD  OF  MISSILE  VELOCITY 

Y-COORD  OF  MISSILE  VELOCITY 

X-COORD  OF  MISSILE  ACCELERATION 

Y-COORD  OF  MISSILE  ACCELERATION 


SHIP  LENGTH 
SHIP  BREADTH 
COURSE 

FIRST  LEG  OF  NEW  COURSE  TO  STEER 
SECOND  LEG  OF  NEW  COURSE  TO  STEER 
COURSE  OPTION  FOR  FIRST/  SECOND  LEG 
DIFFERENCE  BETWEEN  CSE1  AND  CSE2 
TIME  TO  TURN  TO  CSEOPT 
ANGLE  BETWEEN  REF  AND  SHIP'S  HEAD 
SCALED  SHIP  VELOCITY 
SHIP  ACCELERATION 
TURN  RATE  OF  SHIP 
,U3)TURN  DIRECTION  OF  SHIP 

ASPECT  ANGLE  OF  SHIP (BETWEEN  CSE  &  BETA) 
X-COORD  OF  SHIP  POSITION 
Y-COORD  OF  SHIP  POSITION 
X-COORD  OF  SHIP  VELOCITY 
Y-COORD  OF  SHIP  VELOCITY 
X-COORD  OF  SHIP  ACCELERATION 
Y-COORD  OF  SHIP  ACCELERATION 
RADAR  CROSS  SECTION  OF  SHIP 
MAXIMUM  RADAR  CROSS  SECTION  OF  SHIP 
MINIMUM  RADAR  CROSS  SECTION  OF  SHIP 
INITIAL  RADAR  CROSS  SECTION  OF  SHIP 
INSTANTANEOUS  RADAR  CROSS  SECTION  OF  SHIP 
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c 

BETA 

FORTRAN 

& 

C 

C      RESOLUTION  CELL: 

C 

RGATE 

C 

RG 

C 

AG 

C 

RGS 

CENTER 

C 

C 

C 

AGS 

CENTER 

C 

C 

ANGLE  BETWEEN  REF  USED  FOR  MICROSOFT 


C 

CENTER 
C 


S2r/S2a 


BRGM 


RANGE  GATE 
SCALED  RANGE  GATE 
ANGLE  GATE 

DISTANCE  BETWEEN  RESCELL  CENTER  &  SHIP 

ALONG  THE  RANGE  GATE 
RGC 

CHAFF 

DISTANCE  BETWEEN  RESCELL  CENTER  &  SHIP 


c 

Slr/Sla 

c 

c 

PARTS 

c 

PARTC 

c 

CHAFF: 

c 

c 

BRGC 

c 

CDIST 

c 

CRAD 

c 

VC 

c 

GC 

c 

BRGC1&2 

CHAFF 

C 

RCSC 

C 

RCSSC(l) 

C 

RCSSC(N)         I 

C 

XC 

c 

YC 

c 

EQUIVALENT  TARGET: 

c 
c 

BRGTM 

c 

GET 

OR 

C 

C 

XET 

C 

YET 

c 

VET 

TARGET 

c 

ALONG  ANGLE  GATE 
AGC 

CHAFF 

DISTANCE  BETWEEN  RG/AG  BOUNDARY  &  SHIP 

Rlr/Rla 

CHAFF 

PORTION  OF  HALF  OF  THE  LONGITUDINAL  AXIS 
INSIDE  THE  CELL 

PORTION  OF  THE  SHIP  INSIDE  THE  RESCELL 
CHAFF 


BEARING  TO  LAUNCH  CHAFF 

CHAFF  LAUNCH  DISTANCE 

CHAFF  RADIUS 

CHAFF  SPEED 

DIRECTION  OF  MOTION  OF  CHAFF 

BEARINGS  PERPENDICULAR  TO  BRGM  TO  FIRE 

RADAR  CROSS  SECTION  OF  CHAFF 
INITIAL  RADAR  CROSS  SECTION  OF  CHAFF (=0) 
INSTANTANEOUS  RADAR  CROSS  SECTION  OF  CHAFF 
X-COORD  OF  CHAFF  POSITION 
Y-COORD  OF  CHAFF  POSITION 


BEARING  OF  MISSILE  FROM  EQUIVALENT  TARGET 
DIRECTION  OF  MOTION  OF  EQUIVALENT  TARGET 

EQUIVALENT  TARGET  ANGLE 
X-COORD  OF  EQUIVALENT  TARGET  POSITION 
Y-COORD  OF  EQUIVALENT  TARGET  POSITION 
RATE  OF  CHANGE  OR  VELOCITY  OF  EQUIVALENT 

POSITION 


WIND: 


85 


c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


GW  WIND  DIRECTION (FROM  WHERE  IT  COMES  FROM) 

VW  TRUE  WIND  SPEED 

VW1  SCALED  WIND  SPEED 

TW  TRUE  WIND  DIRECTION (TO  WHERE  IT  GOES) 

RELATIVE  COORDINATES 


XR 
YR 

XRP 


YRP 


XRPD(l) 
YRPD(l) 
YRPD(N) 


VELOCITY 

C 

VELOCITY 

C         XRP1 

DIRECTION 

C         YRP1 

DIRECTION 

C         XRP2 

DIRECTION 

C         YRP2 

DIRECTION 


XRPD(N) 


C 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


XRPD1 
YRPD1 
XRPD2 
YRPD2 

XRPDD 

YRPDD 
MISCELLANEOUS: 

SC 

TTG 

ITTG 

XCENT 

YCENT 

NAV 

SIG 


SIGHT (CLOCKWISE) 

C        SIGD 

C 

10     PI  =  4.0  * 
NAV  =3.0 
T=0 

H  =  1.0 
XCENT  =4.1 


X-COORD  OF  RELATIVE  POSITION  OF  SHIP 
Y-COORD  OF  RELATIVE  POSITION  OF  SHIP 
X-COORD  OF  RELATIVE  POSITION  OF  SHIP 
WITH  THE  RELATIVE  COORDINATE 
SYSTEM  ROTATED  TO  MISSILE  BEARING 
Y-COORD  OF  RELATIVE  POSITION  OF  SHIP 
WITH  THE  RELATIVE  COORDINATE 
SYSTEM  ROTATED  TO  MISSILE  BEARING 
INITIAL  X-COORD  OF  RELATIVE  SHIP'S  VELOCITY 
INITIAL  Y-COORD  OF  RELATIVE  SHIP'S  VELOCITY 
INSTANTANEOUS  Y-COORD  OF  REL  SHIP'S 

INSTANTANEOUS  X-COORD  OF  REL  SHIP'S 

DISTANCE  TRAVELLED  BY  SHIP  IN  POS  XRP 

DISTANCE  TRAVELLED  BY  SHIP  IN  POS  YRP 

DISTANCE  TRAVELLED  BY  SHIP  IN  NEG  XRP 

DISTANCE  TRAVELLED  BY  SHIP  IN  NEG  YRP 

RELATIVE  VELOCITY  OF  SHIP  IN  XRP1  DIRECTION 
RELATIVE  VELOCITY  OF  SHIP  IN  YRP1  DIRECTION 
RELATIVE  VELOCITY  OF  SHIP  IN  XRP2  DIRECTION 
RELATIVE  VELOCITY  OF  SHIP  IN  XRP2  DIRECTION 
X-COORD  OF  RELATIVE  ACCELERATION  OF  SHIP 
Y-COORD  OF  RELATIVE  ACCELERATION  OF  SHIP 


SCREEN  SCALE 
TIME  TO  GO 
PROGRAM  RUN  TIME 
X-COORD  OF  SCREEN  CENTER 
Y-COORD  OF  SCREEN  CENTER 
PROPORTIONAL  NAVIGATION  CONSTANT 
ANGLE  BETWEEN  REF  AND 


LINE  OF 


RATE  OF  CHANGE  OF  LOS 


ATAN(l.O) 
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YCENT  =3.0 
SL  =  200.0 
SB  =  20.0 
IFLAGl(l)  = 
IFLAG2(1)  = 
IFLAG3  =  0 


C 

c 


WRITE (*,*)  'INPUT 
READ(*/*)  BRGM(l) 
WRITE ( * , * )  ' INPUT 
READ(*,*)  RANGE 
WRITE ( * , * )  ' INPUT 
READ(*,*)  CSE(l) 
MISSILE  PARAMETERS 


MISSILE  BEARING  (000-359  DEGREES  TRUE)  ' 
INITIAL  MISSILE  RANGE  (<  15000  YARDS)  ' 
SHIP  HEADING  (000-359  DEGREES  TRUE) ' 


WRITE (*,*)  'INPUT  MISSILE  RADAR  BEAMWIDTH  (+/"  1.0  TO  +/-5.0 
♦DEGREES) ' 
READ ( * , * )  THETA 

WRITE(*,*)  'INPUT  THE  RANGE  GATE  (+/-10  TO  +/-100  YARDS)' 
READ(*,*)  RGATE 
RADAR  CROSS  SECTION 


C 

C    

C    SHIP  RCS 

BOWS/QTRS 

WRITE (* 
WRITE (* 
WRITE (* 
WRITE (* 
WRITE (* 
WRITE (* 
WRITE (* 


ASSUMED  MAX  RIGHT  AHEAD/ASTERN  &  ABEAM  AND  MIN  AT 


READ ( * , * ) 
WRITE (* 
WRITE (* 
WRITE (* 
READ ( * , * ) 


C 
C 


WRITE (* 
WRITE (* 
WRITE (* 
READ(*,*)  RCSC 
WIND  PARAMETERS 


*) 
*) 
*) 
*) 
*) 
*) 
*) 


*) 
*) 
*) 


*) 
*) 

*) 


INPUT  RADAR  CROSS  SECTION  OF  SHIP  &  CHAFF' 
FOR  CHAFF  TO  BE  EFFECTIVE  THERE  MUST  BE  AT 
A  1:2  RCS  RATIO  ( CHAFF: SHIP) ' 


LEAST 


INPUT  MAX  RADAR  CROSS 
10000  -  100000  SQUARE 


SECTION 
METERS ' 


OF  SHIP 


RCSSMAX 
INPUT  MINIMUM  RADAR  CROSS  SECTION  OF  SHIP' 
5000  -  50000  SQUARE  METERS' 


RCSSMIN 
INPUT  RADAR  CROSS  SECTION  OF  CHAFF' 
5000  -  20000  SQUARE  METERS  ' 


20 

C 

C 


WRITE (*,*)  'INPUT 
READ(*,*)  GW 

WRITE (*,*)  'INPUT  TRUE  WIND  SPEED 
READ(*,*)  VW 
WRITE (*,*)  'INPUT 
WRITE (*,*)'  RECOMMEND 
READ(*,*)  SC 
CONTINUE 
SCALING 


WIND  DIRECTION  (000-359  DEGREES  TRUE)  ' 

(0.0-50.0  KNOTS)  ' 


SCREEN  SCALE   (YDS:1 
2  00  YARDS  TO  THE 


INCH) 
INCH' 
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C      VS  =  SHIP'S  SCALED  VELOCITY   CONVERTED  TO  YDS/SEC 

VS(1)  =  20./(SC*1.8) 
C      RG  =  SCALED  RANGE  GATE 

RG  =  RGATE/SC 
C      VM  =  SCALED  MISSILE  VELOCITY 

VMM  =  600. 

VM  -  VMM/(SC*1.8) 
C      VW1  -  SCALED  WIND  VELOCITY 

VW1  =  VW/(SC*1.8) 
C      SCALE  THE  RANGE  FOR  PLOTTING 

R(l)  =  RANGE/SC 

RMIN  =  R(l) 
C      TR  IS  THE  TURN  RATE  IN  DEGREES/SEC 

TR  =  5.0 
C      VR  IS  THE  ACCELERATION  IN  KNOTS/SEC*SC 

VR  =  0.5/(SC*1.8) 

RCSS=RCSSMIN+ABS  (  (RCSSMAX-RCSSMIN)  *COS  (2*  (CSE  (1)  -BRGM(l)  )  *PI/180)  ) 

RCSSP(l)  =  RCSS 

RCSCP(l)  =0.0 
C     TTG  =  TIME  TO  GO 
c     

TTG  =  RANGE* 1.8/VMM 

ITTG  =  IFIX(TTG  +10*500./ (VM*SC) ) 
C    MISSILE  ACCELERATION 
c    

AMM=30./SC 

C     SIG   REF.  TO  HORIZONTAL  (EAST) 
c     

SIG  =  (270.  -  BRGM(l) )*PI/180. 
IF  (BRGM(l) .GT.180.)  THEN 

GM  =  BRGM(l)  -  180. 
ELSE 

GM  =  BRGM(l)  +  180. 

END  IF 

C   ******INITIAL  CONDITIONS  ****** 

C    MISSILE  I.  C. 
c    

XM(1)  =  XCENT  +  R(1)*SIN(BRGM(1)*PI/180.) 

YM(1)  =  YCENT  +  R(1)*COS(BRGM(1)*PI/180.) 

XMD(l)  =  VM*SIN(GM*PI/180.) 

YMD(l)  =  VM*COS(GM*PI/180.) 
C    SHIP  I.  C. 
C    

XS(1)  =  XCENT 

YS(1)  =  YCENT 

XSD(l)  =  VS(1)  *  SIN(CSE(1)*PI/180.) 

YSD(l)  =  VS(1)  *  COS(CSE(1)*PI/180.) 

XSDD(l)  =0.0 

YSDD(l)  =0.0 

GS  =  ATAN2(YSD(1) ,XSD(1) ) 
C    FIND  THE  NEW  COURSE  TO  STEER  AND  BEARING  TO  LAUNCH  CHAFF 
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c   

30     WRITE (*,*)'  INPUT  (1)  FOR  AUTOMATIC  CHAFF  LAUNCH   ' 
WRITE (*,*)'  AND  NEW  SHIP  COURSE  ' 

WRITE (*,*)'   ' 

WRITE (*,*)'  INPUT  (2)  FOR  MANUAL  CHAFF  LAUNCH    ' 
WRITE (*,*)'  AND  NEW  SHIP  COURSE  ' 

READ(*,*)  ISHOOT 
IF  ( ISHOOT. EQ.l)  THEN 

CALL  AUTOl  (CSE2,U2) 

U  =  U2 

CALL  AUT02  (CSE3 /U3/CSE2) 

Ul  =  U3 

C    DRAW  THE  INITAL  MISSILE  RESOLUTION  CELL 
c    

WRITE(*,*) '****   MISSILE   DETECTION  ****  ' 

RMT(l)   =  R(l) 

BRGTM(l)  =  BRGM(l) 

XET(l)  =  XS(1) 

YET(l)  =  YS(1) 

CALL  SCREN1 
ELSE  IF  ( ISHOOT. EQ. 2)  THEN 

WRITE(*,*) '****   MISSILE   DETECTION  ****  ' 

RMT(l)   =  R(l) 

BRGTM(l)  =  BRGM(l) 

XET(l)  =  XS(1) 

YET(l)  =  YS(1) 

CALL  SCREN1 

WRITE (*,*) 'INPUT  TRUE  BEARING  TO  LAUNCH  CHAFF' 

READ(*,*)  BRGC 

WRITE (*,*) 'INPUT  FIRST  LEG  OF  NEW  SHIP  COURSE  TO  STEER  ' 

READ(*,*)  CSE2 
40       WRITE(*,*)'  INPUT  (1)  TO  TURN  THE  SHIP  RIGHT  ' 

WRITE(*,*)'  INPUT  (-1)  TO  TURN  THE  SHIP  LEFT  ' 

READ(*,*)  U2 

U  =  U2 

IF  ( (U2.NE.1.0) .AND. (U2.NE.-1.0) )  THEN 
GO  TO  4  0 

END  IF 

WRITE (*,*) 'INPUT  SECOND  LEG  OF  NEW  SHIP  COURSE  TO  STEER  ' 

READ(*,*)  CSE3 
50       WRITE(*,*)'  INPUT  (1)  TO  TURN  THE  SHIP  RIGHT  ' 

WRITE (*,*)'  INPUT  (-1)  TO  TURN  THE  SHIP  LEFT  ' 

READ(*,*)  U3 

Ul  =  U3 

IF  ((U3.NE.1.0) .AND. (U3.NE.-1.0) )  THEN 
GO  TO  50 

END  IF 

ELSE 

GO  TO  30 

END  IF 

C    CHAFF  INITIAL  CONDITIONS 
C    
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C    CDIST  =  DISTANCE  FM  SHIP  CHAFF  IS  LAUNCHED 

C    CRAD  =  CHAFF  RADIUS 

C    CHAFF  IS  FIRED  TO  THE  EDGE  OF  THE  CELL  BUT  THE  DISTANCE  IS 

REDUCED 

C   BY  ONE  CHAFF  RADIUS  TO  CATER  FOR  THE  SHIP'S  MOVEMENT  DURING  THE 

C    BLOOM  TIME  DELAY 

C    VC  =  CHAFF  SPEED  IN  YDS/SEC 

CRAD  =   100. /SC 

CDIST  =  (R(1)*TAN(THETA  *  PI/180))  -  CRAD 

VC  =  VW  /  (SC*1.8) 

GC  =  GW  +  180. 

IF  (GC.GT.360)  THEN 
GC  =  GC  -  360 

ENDIF 

XC(1)  =  XS(1) 

YC(1)  =  YS(1) 
C 
C 
C   ************  MAIN   LOOP  ****************** 

DO  100  N  -  1,ITTG 

IF  (IFLAGl(N-l) .EQ.l)  GOTO  110 
C 
C 
C 

IFLAGl(N)  =  0 
IFLAG2(N)  =  0 
C    RANGES  AND  ANGLES  (WRT  HORIZONTAL  REFERENCE 
c    


R(N)  =  SQRT( (XS(N)-XM(N) ) **2+ (YS (N) -YM(N) )**2) 

RCM  =  SQRT((XC(N)-XM(N) )**2+(YC(N)-YM(N) )**2) 

RMT(N)  =  ( (XM(N)-XET(N) ) **2+ (YM(N) -YET(N) )**2)**0.5 

IF( ( (XET(N)-XM(N) ) .EQ.O.) .AND. ( (YET (N) -YM(N) ) .EQ.O.) ) THEN 

BRGMT(N)  =0.0 
ELSE 

BRGMT(N)  =  ATAN2(YET(N)-YM(N) , XET(N) -XM(N) ) 
END  IF 
IF( ( (XET(N)-XM(N) ) .EQ.O.) .AND. ( (YET(N) -YM(N) ) .EQ.O.) ) THEN 

BRGTM(N)  =0.0 
ELSE 

BRGTM(N)  =  (180./PI)*ATAN2(XM(N)-XET(N) ,  YM  (N) -YET  (N) ) 
END  IF 
IF( ( (XS(N)-XM(N) ) .EQ.O.) .AND. ( (YS (N) -YM(N) ) .EQ.O.) )  THEN 

BRGMS(N)  =0.0 
ELSE 

BRGMS(N)  =  ATAN2(YS(N)-YM(N) ,XS(N)-XM(N) ) 
END  IF 
IF( ( (XC(N)-XM(N) ) .EQ.O.) .AND. ( (YC(N) -YM(N) ) .EQ.O.) )  THEN 

BRGMC  =0.0 
ELSE 

BRGMC  =  ATAN2(YC(N)-YM(N) ,XC(N)-XM(N) ) 
END  IF 
RANGE  AND  ANGLE  GATES 
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AGS  =  R(N)  *  SIN(ABS(BRGMS(N)-BRGMT(N) ) ) 

AGC  =  RCM  *  SIN(ABS(BRGMC-BRGMT(N))) 

RGS  =  ABS(R(N)  *  COS (ABS (BRGMS (N) -BRGMT(N) ) ) -RMT(N) ) 

RGC  =  ABS (RCM  *  COS (ABS (BRGMC-BRGMT(N) ) ) -RMT(N) ) 

AG  =  RMT(N)  *  TAN(THETA*PI/180.) 
C    RCS  MODEL 
C    


RCSS=RCSSMIN+ABS ( (RCSSMAX-RCSSMIN) *COS (2* (CSE (N) -BRGM(N) ) 

*  *PI/180)) 
RCSSP(N)  =  RCSS 

IF  ((RGS.GT. (RG-SL/(SC*2.0))) .OR. 

*  (AGS.GT. (AG-SL/(SC*2.0))))  THEN 

CALL  SHPCLL(N) 

RCSSP(N)  =  RCSS  *  PARTS 
END  IF 
C    CHAFF  RCS  IS  ASSUMED  TO  BE  ZERO  UNTIL  FULL  BLOOM 
c    

IF  (N.GT.5)  THEN 

RCSCP(N)  =  RCSC 

IF  ( (RGC.GT. (RG-CRAD) ) .OR. (AGC.GT. (AG-CRAD) ) )  THEN 

CALL  CHFCLL 

RCSCP(N)  =  RCSC  *  PARTC 

END  IF 

ELSE 

RCSCP(N)  =0.0 

ENDIF 

C    EQUIVALENT  TARGET 
c    


IF  ((RCSSP(N) .EQ.0.0) .AND. (RCSCP(N) .EQ.0.0) )  THEN 

XET(N+1)  =  XET(N) 

YET(N+1)  =  YET(N) 
ELSE 

XET(N+1)  =  (RCSSP(N)*XS(N)+RCSCP(N)*XC(N) )/(RCSSP(N)+ 

*  RCSCP(N)) 

YET(N+1)  =  (RCSSP(N)*YS(N)+RCSCP(N)*YC(N) )/(RCSSP(N)+ 

*  RCSCP(N)) 
END  IF 

C    ANGLE  BETWEEN  HORIZONTAL  REFERENCE  AND  DIRECTION  OF  MOVEMENT 

OF 

C    EQUIVALENT  TARGET 

IF  (N.GE.2)  THEN 

IF(((XET(N)-XET(N-1) ) .EQ.O.) .AND. ( (YET(N) -YET(N-l) ) .EQ.O.)) 

*  THEN 

GET  =0.0 
ELSE 

GET  =  ATAN2(YET(N)-YET(N-1) ,XET (N) -XET (N-l) ) 

END  IF 

C    EQUIVALENT  TARGET  VELOCITY 
C    

VET  =  SQRT( (XET(N)-XET(N-l) ) **2+ ( YET (N) -YET (N-l) )**2)/H 
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ELSE 

GET  =  GS 
VET  =  VS(N) 
END  IF 
C   MISSILE  NAVIGATION  DYNAMICS 
c    

IF(((XET(N)-XM(N)) .EQ.O.) .AND. ( (YET(N) -YM(N) ) .EQ.O.))THEN 

SIG  =0.0 
ELSE 

SIG  =  ATAN2((YET(N)-YM(N)) , (XET(N)-XM(N))) 
END  IF 
IF  (RMT(N) .EQ.O.O)  THEN 

SIGD  =0.0 
ELSE 

SIGD  =  (VET*SIN(GET-SIG)  -VM*SIN(GM-SIG) )/RMT(N) 
END  IF 

AM=SIGD*VM*NAV 
IF(AM.GT.AMM)  THEN 

AM=AMM 
ELSE  IF(AM.LT.-AMM)  THEN 

AM  =  -AMM 

END  IF 

C    TO  TURN  OFF  MISSILE  SEEKER 
c    


IF  (RMT(N) .LT. (250./SC))  THEN 

AM  =  0.0 

IF  (RMT(N-l) .LT.RMT(N) )  THEN 

ITTG  =  N 

IFLAGl(N)  =  1 

END  IF 

END  IF 

C    MISSILE  EQUATIONS 
c    


XMDD(N)=  -AM*SIN(SIG) 
YMDD(N)=  AM*COS(SIG) 
IF( (YMD(N) .EQ.O.) .AND. (XMD(N) .EQ.O.) ) THEN 

GM  =  0.0 
ELSE 

GM  =  ATAN2(YMD(N) ,XMD(N) ) 

END  IF 

XM(N+1)=  XM(N)  +  XMD(N)*H  +  0 . 5*XMDD(N) *H**2 

YM(N+1)=  YM(N)  +  YMD(N)*H  +  0 . 5*YMDD(N) *H**2 

XMD(N+1)=   VM*COS(GM)  +  XMDD(N)*H 

YMD(N+1)=   VM*SIN(GM)  +  YMDD(N)*H 

C    SHIP   EQUATIONS 
c    


IF  (BRGM(N) .LT.0.0)  THEN 

BRGM(N)  =  BRGM(N)  +  360. 
END  IF 

XSDD(N)   =   -VS(N)*(PI*TR/180.)*U1*SIN(CSE(N)*PI/180.)   + 
VR*SIN 

*   (PI*CSE(N)/180.) 
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YSDD(N)      =     VS(N)*(PI*TR/180.)*U1*COS(CSE(N)*PI/180.)  + 

VR*COS 

*  (PI*CSE(N)/180.) 

XSD(N)  =  VS(N)*SIN(PI*CSE(N)/180.) 
YSD(N)  =  VS(N)*COS(PI*CSE(N)/180.) 
IF    (VS(N) .LT. (25./(1.8*SC)))    THEN 

VS(N+1)    =   VS(N)    +   VR*H 
ELSE 

VS(N+1)=25./(1.8*SC) 
END   IF 

XS(N+1)  =  XS(N)  +  XSD(N)*H  +  0.5*XSDD(N) *H**2 
YS(N+1)  =  YS(N)  +  YSD(N)*H  +  0. 5*YSDD(N) *H**2 
IF(((XM(N+1)-XS(N+1)) .EQ.O.) .AND. ( (YM(N+1) -YS(N+1) ) .EQ.O.)) 

*  THEN 

BRGM(N+1)    =0.0 
ELSE 

BRGM(N      +      1) 
(180./PI)*ATAN2(XM(N+1)-XS(N+1) , YM(N+1) -YS (N+l) ) 

END  IF 
C    CHAFF  DYANAMICS  AND  RELATIVE  COORDINATES  EQUATIONS 
c    

IF  (N.LE.5)  THEN 

C    CHAFF'S  ARBITRARY  POSITION  BEFORE  IT  BLOOMS 
c    

XC(N+1)  =  XS(N+1)  +  CDIST*SIN(BRGC*PI/180.) 

YC(N+1)  =  YS(N+1)  +  CDIST*COS(BRGC*PI/180.) 
ELSE 

XC(N+1)  =  XC(N)  +  VC*SIN(GC*PI/180.) 

YC(N+1)  =  YC(N)  +  VC*COS(GC*PI/180.) 

ENDIF 

C    SHIFT  OF  THE  ORIGIN  OF  COORDINATE  SYSTEM  TO  THE  CHAFF 
c    

XR(N)  =  XS(N)-XC(N) 

YR(N)  =  YS(N)-YC(N) 
C    ROTATION  OF  THE  COORDINATE  SYSTEM  TO  THE  MISSILE  DIRECTION 
C    

XRP(N)  =  (XR(N)*COS(BRGM(N)*PI/180) )  -  (YR(N)* 

*  SIN(BRGM(N)*PI/180) ) 

YRP(N)  =  (XR(N)*SIN(BRGM(N)*PI/180) )  +  (YR(N)* 

*  COS(BRGM(N)*PI/180) ) 
T=T+H 

C    ALTERATION  OF  COURSES 
C    

C    WATCH  FOR  CSE(N)  -  CSE2/CSE3  NEVER  BEING  LESS  THAN  TR 
C    WHILE  THE  SHIP  IS  WITHIN  THE  RG  BOUNDS  UPDATE  COURSE  TO  ALTER 
TO 

C    THE  FIRST  LEG  UNTIL  THE  TURN  IS  COMPLETED 
IF  (ABS(YRP(N) ) .LE.RG)  THEN 

IF  ( (ABS(CSE(N)-CSE2) ) .LT.TR)  THEN 

U  =  0.0 
END  IF 
CSE(N+1)  =  CSE(N)  +  TR*H*U 
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C    WHILE  THE  SHIP  IS  WITHIN  THE  RG  BOUNDS  BUT  IT  HAS  EXITED  THE 

AG 

C    BOUNDS  UPDATE  COURSE  TO  ALTER  TO  THE  FIRST  LEG  UNTIL  THE  TURN 

C    IS  COMPLETED 

IF  (ABS(XRP(N)) .GT.AG)  THEN 

IF  ((ABS(CSE(N)-CSE3)) .LT.TR)  THEN 
Ul  =  0.0 
ENDIF 

CSE(N+1)  =  CSE(N)  +  TR*H*U1 
ENDIF 
ELSE 
C    IF  THE  SHIP  IS  OUT  OF  THE  RG  BOUNDS  UPDATE  COURSE  TO  ALTER  TO 
THE 
C    SECOND  LEG  UNTIL  THE  TURN  IS  COMPLETED 

IF  ((ABS(CSE(N)-CSE3)) .LT.TR)  THEN 
Ul  -  0.0 
ENDIF 

CSE(N+1)  =  CSE(N)  +  TR*H*U1 
ENDIF 
IF  (CSE(N+1) .GE.360.)  THEN 

CSE(N+1)  =  CSE(N+1)  -  360.0 
END  IF 
IF  (CSE(N+1) .LT.O.)  THEN 

CSE(N+1)  =  360.+CSE(N+l) 
END  IF 
C    DETERMINE  IF  MISSILE  HITS  THE  SHIP 
c    

IF( (R(N)*SC) .LT. 50.0) THEN 

WRITE (*,*)  'RANGE  UNDER  50  YARDS  MISSILE  HITS  SHIP' 
END  IF 

RMIN  =  AMIN1(RMIN,R(N) ) 
IF  (RMT(N) .LT. (500./SC) )  THEN 

IFLAG2(N)  =  1 

H  =  0.1 
END  IF 
C 
C 

c 

100   CONTINUE 

110   CONTINUE 

C    TO  PREVENT  ANY  UNDEFINED  VALUES  IF  ITTG  <  60 
c    


DO  120  I=(ITTG+1) ,60 
XC(I)  =  XC(I-l) 
XC(I)  =  XC(I-l) 
YC(I)  =  YC(I-l) 
XS(I)  =  XS(I-l) 
YS(I)  =  YS(I-l) 
XM(I)  =  XM(I-l) 
YM(I)  =  YM(I-l) 
XET(I)  =  XET(I-l) 
YET(I)  =  YET(I-l) 
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120   CONTINUE 
C 
C 
C 

RMIN  =  SC*RMIN 
C* ********************************************************** * 

C    PRINT  &  PLOT 


c 

130 

WRITE(*,*) 'INPUT  (1)  TO  DRAW  MISSILE  &  SHIP  &  CHAFF  PATH' 

WRITE (*,*)'     ' 

WRITE (*,*) 'INPUT  (2)  TO  DRAW  RESOLUTION  CELL  IN  ' 

WRITE (*,*)'           FIVE  SECOND  INCREMENTS      ' 

WRITE (*,*)'     ' 

WRITE (*,*)' INPUT  ( 3 )  TO  CHANGE  SCENARIO  PARAMETERS ' 

WRITE (*,*)'     ' 

WRITE (*,*) 'INPUT  (0)  TO  QUIT         ' 

WRITE (*,*)'     ' 

READ(*,*)  IDRW 

IF(IDRW.EQ.l)  THEN 

CALL  PATH 

ELSE  IF(IDRW.EQ.2)  THEN 

CALL  SCREN2 

ELSE  IF(IDRW.EQ.3)  THEN 

CALL  SCREEN3(IFLAG3) 

IF  (IFLAG3.EQ.1)  THEN 

GO  TO  10 

ENDIF 

IF  (IFLAG3.EQ.2)  THEN 

PI  =  4.0  *  ATAN(l.O) 

NAV  =3.0 

T=0 

H  =  1.0 

XCENT  =4.0 

YCENT  =3.0 

SL  =  200.0 

SB  =  20.0 

IFLAGl(l)  =  0 

IFLAG2(1)  =  0 

GO  TO  20 

ENDIF 

ELSE  IF(IDRW.EQ.O)  THEN 

GO  TO  999 

END  IF 

GO  TO  13  0 

999 

STOP 

C 

END 

c**************************************************************** 

SUBROUTINE  AUTOl  (CSE2,U2) 

c 

THIS  SUBROUTINE  CALCULATES  THE  BEARING  TO  LAUNCH  THE  CHAFF  AND 

THE 

C 

FIRST  LEG  OF  THE  NEW  COURSES  TO  MANEUVER  THE  SHIP 
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COMMON/TIME/ ITTG , TTG 
COMMON/WD/ GW , VW 

COMMON/ABCD/BRGM(60) ,R(60) ,SC,VS (60) ,CSE(60) ,BRGMT(60) , 
*BRGTM(60) ,RMT(60) 
COMMON/MISS/RG , THETA , AG , PI 
COMMON/ CHF/BRGC, I SHOOT 

COMMON/REL/XR(60)  ,YR(60)  ,XRD(60)  ,YRD(60)  ,XRP(60)  ,YRP(60)  ,XRPD(60)  , 

*YRPD(60) 

COMMON/TURN/TR, CSEOFF 

REAL  CSEOPT ( 4 ) , TOPT (  4 )  , CSE2 , CSEDIF ( 4 ) , CDIF1 , ACC , VSMAX 

REALXRP1 , YRP1 , XRP2 , YRP2 , XRPD1 , YRPD1 , XRPD2 , YRPD2 , XRPDD, YRPDD 

INTEGER  TTG1,T 

REAL  CDIST 

C    POSSIBLE  CHAFF  LAUNCH  BEARINGS 
c    

BRGC1  =  BRGM(l)  +  90 
BRGC2  =  BRGM(l)  -  90 
IF  (BRGC1.GT.360)  THEN 

BRGC1  =  BRGC1  -  360 
ELSEIF  (BRGC1.LT.0)  THEN 

BRGC1  =  BRGC1  +  360 
ENDIF 
IF  (BRGC2.GT.360)  THEN 

BRGC2  =  BRGC2  -  360 
ELSEIF  (BRGC2.LT.0)  THEN 

BRGC2  =  BRGC2  +  3  60 

ENDIF 

WRITE  (*,*)     'BRGC1=',BRGC1 

WRITE  (*,*)  'BRGC2=',BRGC2 

GC  =  GW  +  180 

VW1  =  VW/(SC*1.8) 

C    COURSE  OPTIONS  FOR  THE  FIRST  LEG 
c    

CSEOPT (1)  =  BRGM(l) 

CSEOPT (2)  =  BRGM(l)  +  180 

C    COURSE  OPTIONS  FOR  THE  SECOND  LEG 
c    

CSEOPT (3)  =  BRGM(l)  +  90 

CSEOPT (4)  =  BRGM(l)  -  90 

DO  10  1=1,4 

IF  (CSEOPT(I) .GT.360)  THEN 

CSEOPT(I)  =  CSEOPT(I)  -  360 

ELSEIF  (CSEOPT (I) .LT.O)  THEN 

CSEOPT (I)  =  CSEOPT (I)  +360 

ENDIF 

10     CONTINUE 

C    CALCULATION  OF  INITIAL  VELOCITY  OF  RELATIVE  VELOCITY 
c    

TEMP1=(VS(1)*(SIN( (CSE(l)-BRGM(l) )*PI/180) ) ) 
TEMP2=(VW1*(SIN( (GW-BRGM(l) )*PI/180) ) ) 
XRPD(l)  =  TEMP1+TEMP2 
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TEMP3=(VS(1)*(COS((CSE(1)-BRGM(1))*PI/180))) 

TEMP4= (VW1* (COS ( (GW-BRGM(l) ) *PI/180) ) ) 

YRPD(l)  =  TEMP3+TEMP4 
C   COURSE  DIFFERENCE  BETWEEN  INITIAL  COURSE  AND  THE  COURSE  OPTIONS 
c 

DO  20  I  =  1,4 

CSEDIF(I)  =  ABS(CSE(1)  -  CSEOPT(I)) 
IF  (CSEDIF(I) .GT.180)  THEN 

CSEDIF(I)  =  360  -  CSEDIF(I) 
ENDIF 
20     CONTINUE 

C    TIME  TAKEN  TO  TURN  TO  COURSE  OPTIONS  3  &  4 
c    _ 


TOPT(3)  =  CSEDIF(3)/TR 

TOPT(4)  =  CSEDIF(4)/TR 
C    INITIAL  VALUES  OF  VARIABLES  INVOLVED 
c    

VSMAX  =  25/(1.8*SC) 

TTG1  =  IFIX(TTG) 

ACC  =  VSMAX  *  TR  *  PI/ 180 

T  =  1 

CRAD  =  100. /SC 

CDIST  =  (R(l)  *  TAN(THETA*PI/180) )  -  CRAD 
C    PRE-COMPUTATION  TO  CALCULATE  CHAFF  LAUNCH  BEARING 
C    

XRP1  =  CDIST 
YRP1  =0.0 
XRP2  =  -XRP1 
YRP2  =0.0 


XRPD1  =  XRPD(l) 

YRPD1  =  YRPD(l) 

XRPD2  =  XRPD(l) 

YRPD2  =  YRPD(l) 

XRPDD  =0.0 

YRPDD  =0.0 

C 

IF  COURSE  IS  ALTERED  TOWARDS  COURSE  OPTION 

C 

3 (TO  THE  RIGHT) 


DO  30  I  =  1,TTG1 

IF  (I.LT.INT(TOPT(3) ) )  THEN 

XRPDD  =  ACC  *  COS(ABS(CSE(l)  -  BRGM(l))  *  PI/180) 

YRPDD  =  -ACC  *  SIN(ABS(CSE(1)  -  BRGM(l) )  *  PI/180) 
ELSE 

XRPDD  =0.0 

YRPDD  =0.0 

ENDIF 

XRP1  =  XRP1  +  XRPD1*T  +  0.5  *  XRPDD*T**2 

YRP1  =  YRP1  +  YRPD1*T  +  0.5  *  YRPDD*T**2 

XRPD1  =  XRPD1  +  XRPDD 

YRPD1  =  YRPD1  +  YRPDD 

30    CONTINUE 

C    IF  COURSE  IS  ALTERED  TOWARDS  COURSE  OPTION  4  (TO  THE  LEFT) 
C    

97 


c 
c 


DO  40  I=1,TTG1 

IF  (I.LT.INT(TOPT(4)))  THEN 

XRPDD  =  -ACC  *  COS(ABS(CSE(l)  -  BRGM(l) )  *  PI/180) 
YRPDD  =  ACC  *  SIN(ABS(CSE(1)  -  BRGM(l))  *  PI/180) 


ELSE 

XRPDD  =0.0 

YRPDD  =0.0 

ENDIF 

XRP2  =  XRP2  +  XRPD2*T  + 

0.5  *  XRPDD*T**2 

YRP2  =  YRP2  +  YRPD2*T  + 

0.5  *  YRPDD*T**2 

XRPD2  =  XRPD2  +  XRPDD 

YRPD2  =  YRPD2  +  YRPDD 

40 

CONTINUE 

C 
C 

CALCULTAION  OF  CHAFF  LAUNCH 

BEARING 

IF  (ABS(XRPl) .GT.ABS(XRP2))  THEN 

BRGC  =  BRGC2 
ELSE 

BRGC  ■  BRGC1 
ENDIF 
CALCULATION  OF  THE  FIRST  LEG  OF  THE  NEW  COURSES 


IF  (YRPD(l) .GT.O)  THEN 
CSE2  =  CSEOPT(l) 

ELSE 

CSE2  =  CSEOPT(2) 

ENDIF 

WRITE  (*,*) 'CSE2=',CSE2 
COURSE  DIFFERENCE  BETWEEN 


INITIAL  COURSE  AND  THE  FIRST  LEG 


CDIF1  =  ABS(CSE(1)  -  CSE2) 

C    CALCULATION  OF  TURN  DIRECTION  INDICATOR  FOR  THE  FIRST  LEG 
c    

IF  (CDIF1.GE.180)  THEN 

IF  (CSE2.GT.CSE(1) )  THEN 

U2  =  -1.0 
ELSE 

U2  =  1.0 
ENDIF 
ELSEIF  (CDIF1.LT.180)  THEN 
IF  (CSE2.GT.CSE(1) )  THEN 

U2  =  1.0 
ELSE 

U2  =  -1.0 
ENDIF 
ENDIF 
RETURN 
END 
C  *************************************************************** 

SUBROUTINE  AUT02  (CSE3 ,U3 , CSE2) 
C    THIS  SUBROUTINE  CALCULATES  THE  SECOND  LEG  OF  THE  NEW  COURSES 
C    TO  MANEUVER  THE  SHIP 
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COMMON/MI SS/RG , THETA , AG , PI 
COMMON/ CHF/BRGC, I SHOOT 

COMMON/REL/XR(60)  ,YR(60)  ,XRD(60)  ,YRD(60)  ,XRP(60)  ,YRP(60)  ,XRDP(60)  , 

*YRDP(60) 

COMMON/ABCD/BRGM(60) ,R(60) ,SC, VS (60) ,CSE(60) ,BRGMT(60) , 

*BRGTM(60) ,RMT(60) 

REAL  CSEOPT ( 4 ) , CSE3 , BRGM 

CSEOPT(l)  =  BRGM(l) 

CSEOPT(2)  =  BRGM(l)  +  180 

CSEOPT (3)  =  BRGM(l)  +  90 

CSEOPT(4)  =  BRGM(l)  -  90 

DO  10  1=1,4 

IF  (CSEOPT(I) .GT.360)  THEN 

CSEOPT(I)  =  CSEOPT(I)  -  360 

ELSEIF  (CSEOPT (I) .LT.O)  THEN 

CSEOPT(I)  =  CSEOPT(I)  +  360 

ENDIF 

10     CONTINUE 

BRGC1  =  BRGM(l)  +  90 

BRGC2  =  BRGM(l)  -  90 

IF  (BRGC1. GT.3  60)  THEN 

BRGC1  =  BRGC1  -  360 

ELSEIF  (BRGC1.LT.0)  THEN 

BRGC1  =  BRGC1  +3  60 

ENDIF 

IF  (BRGC2. GT.360)  THEN 

BRGC2  =  BRGC2  -  360 

ELSEIF  (BRGC2.LT. 0)  THEN 

BRGC2  =  BRGC2  +  3  60 

ENDIF 

C    CALCULATION  OF  THE  SECOND  LEG  OF  THE  NEW  COURSES 
c    

IF  (BRGC.EQ.BRGC1)  THEN 

CSE3  =  CSEOPT (4) 

ELSE 

CSE3  =  CSEOPT (3) 

ENDIF 

WRITE (*,*) 'CSE3=',CSE3 

C    CALCULATION  OF  THE  TURN  DIRECTION  INDICATOR  FOR 

C    THE  SECOND  LEG  OF  THE  NEW  COURSES 

C    

IF  ( CSE3.EQ. CSEOPT (3 ) )  THEN 
IF  ( CSE2.EQ. CSEOPT ( 1) )  THEN 

U3  =  1.0 
ELSEIF  ( CSE2.EQ. CSEOPT (2 ) )  THEN 

U3  =  -1.0 
ENDIF 
ELSEIF  (CSE3.EQ.CSEOPT(4) )  THEN 
IF  ( CSE2.EQ. CSEOPT ( 1) )  THEN 

U3  =  -1.0 
ELSEIF  (CSE2. EQ. CSEOPT (2 ) )  THEN 
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U3  =  1.0 
ENDIF 
ENDIF 
RETURN 
END 
C* *************************************************************** 

******** 

SUBROUTINE  PATH 
C    THIS  SUBROUTINE  OF  PLOT-88  DISPLAYS  THE  PATHS  FOLLOWED 
C    BY  THE  SHIP,  MISSILE  AND  CHAFF 

COMMON/ PRN/IOPORT , MODEL 

COMMON/TIME/ ITTG , TTG 

COMMON/TRACK/  XM(60) ,YM(60) ,XS(60) ,YS(60) ,XC(60) ,YC(60) 

COMMON/MINN/RMIN 

REAL  XCP(55) ,YCP(55) 

CHARACTER*20  LS,LM,LC, YD 

CHARACTER* 2 5  LRMIN 

YD  =  'YARDS' 

LRMIN= 'MINIMUM  RANGE  TO  SHIP' 

C    SHIFT  OF  CHAFF  POSITION  ARRAY  TO  PLOT  CHAFF  AFTER  IT  BLOOMS 
c    

DO  10  I  =  1, (ITTG-5) 

XCP(I)  =  XC(I+5) 

YCP(I)  =  YC(I+5) 

10     CONTINUE 

C    TO  PREVENT  UNDEFINED  VALUES  IF  ITTG-5  <  55 
c    

DO  20  I=(ITTG-4) ,55 
XCP(I)  =  XCP(I-l) 
YCP(I)  =  YCP(I-l) 
20    CONTINUE 

LS  =  'S' 

LM  =  'M' 

LC  =  'C 

K  =  2*ITTG/3 

CALL  PLOTS (0,IOPORT, MODEL) 

CALL  CURVE (XM,YM, -ITTG, 0.1) 

CALL  SYMBOL (XM(K) , YM (K),0. 20, LM, 0.0,1) 

CALL  CURVE (XS,YS, -ITTG, 0.1) 

CALL  SYMBOL(XS(l) , YS (1),0. 20, LS, 0.0,1) 

CALL  CURVE (XCP,YCP, -ITTG, 0.1) 

CALL  SYMBOL(XC(6) ,YC(6) , 0. 20, LC, 0. 0, 1) 

CALL  SYMBOL( 0.5, 0.2, 0.25, LRMIN, 0.0, 21) 

CALL  NUMBER (5.5,0.2,0.25, RMIN ,0.0,1) 

CALL  SYMBOL (7. 0,0. 2, 0.25, YD, 0.0, 5) 

CALL  PLOT(0, 0,999) 

RETURN 

END 
C**************************************************************** 

SUBROUTINE  SCREN1 
C    THIS  PLOT-88  SUBROUTINE  DISPLAYS  THE  INITIAL  CONDITIONS 
C    OF  THE  SHIP  AND  THE  MISSILE 
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COMMON/ PRN/IOPORT , MODEL 

COMMON/ CHF/BRGC, ISHOOT 

CHARACTER* 2 5  CL 

CL  =  'AUTO  CHAFF  LAUNCH  BEARING' 

CALL  PLOTS ( 0 , IOPORT , MODEL) 

CALL  SHIP(l) 

CALL  MISSLE(l) 

CALL  RESCLL(l) 

CALL  WIND 

CALL  TABLEl(l) 

IF  ( ISHOOT. EQ.l)  THEN 

CALL  SYMBOL(l. 0,5.7, 0. 16,CL, 0. 0,25) 

CALL  NUMBER (6.0,5.7,0.16, BRGC, 0 .0,-1) 
END  IF 

CALL  PLOT(0, 0,999) 
RETURN 
END 
C********************************************* ******************* 

SUBROUTINE  SCREN2 
C    THIS  PLOT-88  SUBROUTINE  DISPLAYS  THE  INSTANTANEOUS  POSITIONS 
C    OF  THE  SHIP,  CHAFF  AND  THE  MISSILE  EVERY  5  SECONDS 
COMMON/ PRN/IOPORT , MODEL 
COMMON/ TIME/ I TTG , TTG 
COMMON/FG/IFLAG1(60) ,IFLAG2(60) 
L  =  ITTG/5 
DO  31  J=1,L 
K  =  J  *  5 
IF  (IFLAG2(N) .EQ.l)  THEN 

K  =  J+l 
END  IF 

CALL  PLOTS (0, IOPORT, MODEL) 
CALL  SHIP(K) 
CALL  MISSLE(K) 
IF  (IFLAGl(K) .EQ.O)  THEN 

CALL  RESCLL(K) 
END  IF 
CALL  WIND 
C    CHAFF  NOT  DISPLAYED  FOR  THE  FIRST  5  SECONDS  TO  SHOW  BLOOM  TIME 
DELAY 
C 


IF  (J.GT.l)  THEN 

CALL  CHAFF (K) 
ENDIF 

CALL  TABLEl(K) 
CALL  PLOT (0,0,999) 
31     CONTINUE 
RETURN 
END 
C**************************************************************** 

SUBROUTINE  SHIP (J) 
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c 
c 


c 
c 


THIS  SUBROUTINE  CALCULATES  THE  COORDINATES  OF  THE  POLYGON 
REPRESENTING  THE  SHIP 

COMMON/TRACK/  XM(60) ,YM(60) ,XS(60) ,YS(60) ,XC(60) ,YC(60) 

COMMON/MI SS/RG , THETA , AG , PI 

COMMON/ABCD/BRGM(60) ,R(60) ,80,78(60) ,CSE(60) ,BRGMT(60) , 
*BRGTM(60) /RMT(60) 

COMMON/SHP/SL, SB 

REAL  XSHIP(5) ,YSHIP(5) 
ROTATES  THE  SHIP  DEPENDENT  UPON  THE  COURSE 


C 
C 


XSHIP(1)=XS 

XSHIP(2)=XS 
♦(SB/2.) ♦SIN 

XSHIP(3)=XS 
*(SB/2.)*SIN 

XSHIP(4)=XS 
*(SB/2.)*SIN 

XSHIP(5)=XS 
*(SB/2.)*SIN 

YSHIP(1)=YS 

YSHIP(2)=YS 
♦(SB/2.) ♦COS 

YSHIP(3)=YS 
♦(SB/2.) ♦COS 

YSHIP(4)=YS 
♦(SB/2.) ♦COS 

YSHIP(5)=YS 

*(SB/2.)*COS 

THIS  PLOT-88 


J)  +  (SL/2. 
J)  +  (SL/4. 
(CSE(J)+90. 
J)  -  (SL/2. 
(CSE(J)+90. 
J)  -  (SL/2. 
(CSE(J)-90. 
J)  +  (SL/4. 
(CSE(J)-90. 
J)  +  (SL/2. 
J)  +  (SL/4. 
(CSE(J)+90. 
J)  -  (SL/2. 
(CSE(J)+90. 
J)  -  (SL/2. 
(CSE(J)-90. 
J)  +  (SL/4. 
(CSE(J)-90. 


*SIN(CSE(J)*PI/180.)/SC 

*SIN (CSE (J) *PI/180 . ) /SC 

*PI/180.)/SC 

♦SIN(CSE(J) ♦PI/180. )/SC 

♦PI/180. )/SC 

*SIN(CSE(J) *PI/180. )/SC 

♦PI/180. )/SC 

*SIN(CSE(J)*PI/180.)/SC 

♦PI/180. )/SC 

♦COS (CSE(J) ♦PI/180. )/SC 

♦COS (CSE(J) ♦PI/180. )/SC 

♦PI/180. )/SC 

♦COS (CSE(J) ♦PI/180. )/SC 

♦PI/180. )/SC 

♦COS (CSE(J) ♦PI/180. )/SC 

♦PI/180. )/SC 

♦COS(CSE(J) ♦PI/180. )/SC 

♦PI/180. )/SC 


SUBROUTINE  FILLS  THE  SHIP  POLYGON 


CALL  STFILL(l) 

CALL  FILL  (xship,yship,5) 

RETURN 

END 

SUBROUTINE  RESCLL(J) 
C    THIS  SUBROUTINE  CALCULATES  THE  INSTANTANEOUS  COORDINATES  OF 
C    THE  RADAR  RESOLUTION  CELL 

COMMON/MI SS/RG , THETA , AG , PI 

COMMON/ABCD/BRGM(60) , R ( 60) , SC, VS (60) , CSE (60) ,BRGMT(60) , 
♦BRGTM(60) ,RMT(60) 

COMMON/EQ/  XET(60) ,YET(60) 

REAL  XCELL(4) ,YCELL(4) 
C    SET  UP  THE  RESOLUTION  CELL  OF  THE  MISSILE 


RMBD  =  BRGTM(J)  +  18  0. 

AG  =  RMT(J)  ♦  TAN(THETA^PI/180.) 

ALPHA  =  90.  +  RMBD 

XW1  =  XET(J)  +  AG*SIN(ALPHA^PI/180.) 

YW1  =  YET (J)  +  AG*COS(ALPHA*PI/180.) 

XW2  =  XET(J)  -  AG*SIN(ALPHA*PI/180.) 

YW2  =  YET (J)  -  AG*COS(ALPHA^PI/180.) 
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c 
c 


XCELL(l)  =  XW1  + 

YCELL(l)  =  YW1  + 

XCELL(2)  =  XW1  - 

YCELL(2)  =  YW1  - 

XCELL(3)  =  XW2  - 

YCELL(3)  =  YW2  - 

XCELL(4)  =  XW2  + 

YCELL(4)  =  YW2  + 
DRAW  THE  RESOLUTION 


RG*SIN (RMBD*PI/180 . 
RG*COS (RMBD*PI/180 . 
RG*SIN (RMBD*PI/180 . 
RG*COS (RMBD*PI/180 . 
RG*SIN (RMBD*PI/180 . 
RG*COS (RMBD*PI/180 . 
RG*SIN (RMBD*PI/180 . 
RG*COS (RMBD*PI/180 . 
CELL 


CALL  STFILL(3) 
CALL  FILL  (XCELL,YCELL,4) 
RETURN 
END 
C* *************************************************************** 

SUBROUTINE  WIND 
C    THIS  PLOT-88  SUBROUTINE  DISPLAYS  THE  WIND  DIRECTION  LINE 

COMMON/MISS/RG , THETA , AG , PI 

COMMON/ WD/ GW , VW 

CHARACTER* 20  LW 

LW  ='W' 

XCENT  =4.0 

YCENT  =3.0 
C      DRAW  WIND  DIRECTION  LINE 

TWBRG  =  GW  +  180. 

XWDP1  =  XCENT  -  2.8*SIN(TWBRG*PI/180.) 

YWDP1  =  YCENT  -  2 . 8*COS (TWBRG*PI/180 . ) 

XWDP2  =  XCENT  -  2 . 0*SIN (TWBRG*PI/180 . ) 

YWDP2  =  YCENT  -  2 . 0*COS (TWBRG*PI/180 . ) 

CALL  PLOT (XWDP1,YWDP 1,3) 

CALL  PLOTD ( XWDP2 , YWDP2 , 2 ) 

CALL  SYMBOL ( XWDP1 , YWDP1 ,0. 20, LW, 0.0,1) 

RETURN 

END 
C**************************************************************** 

SUBROUTINE  CHAFF (J) 
C    THIS  SUBROUTINE  CALCULATES  THE  COORDINATES  THE  COORDINATES 
C    OF  THE  HEXAGON  REPRESENTING  THE  CHAFF 

COMMON/TRACK/  XM(60) ,YM(60) ,XS(60) ,YS(60) ,XC(60) ,YC(60) 

COMMON/ABCD/BRGM(60) , R ( 60) , SC, VS (60) ,CSE(60) ,BRGMT(60) , 
*BRGTM(60) ,RMT(60) 

COMMON/ CFCL/ PARTC , CRAD , RGC , AGC 

COMMON/MISS/RG , THETA , AG , PI 

REAL  XCF(6) ,YCF(6) 


c 
c 

DRAW  CHAFF  CLOUD 

CRAD 

XCF(1)=XC(J)  + 

YCF(1)=YC(J) 

XCF(2)=XC(J)  + 

CRAD 

* 

YCF(2)=YC(J)  + 

CRAD 

* 

XCF(3)=XC(J)  - 

CRAD 

* 

COS(60.*PI/180.) 
SIN(60.*PI/180.) 
COS(60.*PI/180. ) 
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YCF(3)=YC(J)  +  CRAD  *  SIN (60. *PI/180. ) 
XCF(4)=XC(J)  -  CRAD 
YCF(4)=YC(J) 

XCF(5)=XC(J)  -  CRAD  *  COS (60. *PI/180. ) 
YCF(5)=YC(J)  -  CRAD  *  SIN(60. *PI/180. ) 
XCF(6)=XC(J)  +  CRAD  *  COS (60. *PI/180. ) 
YCF(6)=YC(J)  -  CRAD  *  SIN(60. *PI/180. ) 
CALL  STFILL(2) 
CALL  FILL  (XCF,YCF,6) 
RETURN 
END 
C* *************************************************************** 

SUBROUTINE  MISSLE(J) 
C    THIS  PLOT-88  SUBROUTINE  DISPLAYS  THE  MISSILE  BEARING  /LINE  OF 
FIRE 

COMMON/EQ/  XET(60) ,YET(60) 

COMMON/MISS/RG , THETA , AG , PI 

COMMON/TRACK/  XM(60) ,YM(60) ,XS(60) ,YS(60) ,XC(60) ,YC(60) 

COMMON/ABCD/BRGM(60) ,R(60) ,SC,VS(60) /CSE(60) ,BRGMT(60) , 
*BRGTM(60) ,RMT(60) 

COMMON/FG/IFLAG1(60) , IFLAG2(60) 

CHARACTER* 20  LM 

LM  =  'M' 

XM1  =  XET(J)  +  2.5  *  SIN(BRGTM(J)*PI/180.) 

YM1  =  YET (J)  +  2.5  *  COS (BRGTM(J) *PI/180. ) 

C    DRAW  THE  MISSILE  BEARING  LINE 
c    

IF  (IFLAGl(J) .EQ.O)  THEN 
IF  (RMT(J) .GT.2.5)  THEN 

CALL  PLOT(XET(J) ,YET(J) ,3) 
CALL  PL0T(XM1,YM1,2) 
CALL  SYMB0L(XM1, YM1 , 0 . 20, LM, 0.0,1) 
ELSE 

CALL  PLOT(XM(J) ,YM(J) ,3) 
CALL  PLOT(XET(J) ,YET(J) ,2) 
CALL  SYMBOL(XM(J) ,YM(J) , 0.20,LM, 0. 0, 1) 
END  IF 
ELSE 

CALL  SYMBOL (XM (J) ,YM(J),0.20,LM,0.0,1) 
END  IF 
RETURN 
END 
C**************************************************************** 

SUBROUTINE  TABLEl(J) 
C     THIS  PLOT-88  SUBROUTINE  DISPLAYS  THE  VALUES  OF  IMPORTANT 
VARIABLES 

C    BOTH  FOR  THE  INITIAL  PICTURE  AND  FOR  THE  RUN  OF  PLAY  EVERY  5 
SECS 

COMMON/ TME/T 

COMMON/ABCD/BRGM(60) , R (60) , SC, VS (60) ,CSE(60) ,BRGMT(60) , 
*BRGTM(60) ,RMT(60) 
COMMON/WD/ GW,VW 
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COMMON/ SHCL/ AGS , RGS , BRGMS ( 60) , SAA, PARTS , SAAN ( 60 ) 
COMMON/RC/RCSS,RCSSP(60)  ,RCSC,RCSCP(60) 

CHARACTER 
CS , SPD , WND , MLL , RGE , SCL , TR , AA , SRCS , CRCS , SP , BR , SCR , TM 


TM   = 

'#' 

BR   = 

'BRG' 

SP   = 

'SHIP' 

CS   = 

'CSE' 

SPD  = 

'SPD' 

MLL  = 

'MISSILE' 

RGE  = 

'RG' 

WND  = 

'WIND' 

SCR  = 

'SCREEN' 

SCL  = 

'SCALE' 

TR   = 

'TRUE' 

SRCS  = 

=  'SHIP   RCS' 

CRCS  = 

-■    'CHAFF  RCS' 

c 
c 

CONVERT 

BACK  TO  THE  UNSCALED  RANGE  &  SPEED 

ZJ    =   J 

RR   =   R(J)     *    SC 

SHSPD   =    1.8    *    SC    *    VS(J) 


c 
c 

DRAW  TABLE  AT  RIGHT  OF  THE  SCREEN 

CALL  PLOT (6. 9 ,6. 0,3) 

CALL  PLOT (6. 9, 0.0, 2) 

call 

symbol (7.0,5.8, 

,0.18,SP,0.0,4) 

call 

symbol (7 . 0 

,5.5 

,  0.16,CS,0.0,3) 

call 

number (7 . 6 

,5.5, 

0.16,CSE(J) ,0.0,-1) 

call 

symbol (7 . 0 

,5.3, 

,0.16, SPD, 0.0, 3) 

call 

number (7 . 6 

,5.3, 

0.16, SHSPD, 0.0,-1) 

call 

symbol (7.0 

,4.8 

0.16,WND,0.0,4) 

call 

symbol (7.0 

,4.5, 

0.16, BR, 0.0, 3) 

call 

number (7 . 6 

,4.5 

,0. 16, GW, 0.0,-1) 

call 

symbol (7.0 

,4.0, 

,0.16, SPD, 0.0, 3) 

call 

number (7 . 6 

,4.0, 

,0. 16, VW, 0.0,-1) 

call 

symbol (7.0 

,3.5, 

,0. 16, MLL, 0.0, 7) 

call 

symbol (7.0 

,3.0, 

,0.16, BR, 0.0, 7) 

call 

number (7. 6 

,3.0, 

0.16,BRGM(J) ,0.0,-1) 

call 

symbol (7.0 

,2.5, 

,0. 16, RGE, 0.0, 2) 

call 

number (7 . 4 

,2.5, 

0. 16, RR, 0.0,-1) 

call 

symbol (7. 0 

,1.7 

0.16, SCR, 0.0, 6) 

call 

symbol (7.2 

,1.4, 

,0.16, SCL, 0.0, 5) 

call 

number (7 . 3 

,1.1, 

0.16,SC, 0.0,1) 

call 

symbol (7.1 

,0.8 

,0.16,TM,0.0,3) 

call 

number (7 .4 

,0.8 

,0. 16, ZJ, 0.0,-1) 

c 
c 

SHIP  RCS 

call 

symbol (4.0,0.1 

,0. 16, SRCS, 0.0, 9) 

call 

number (6. 0,0.1 

,0.16,RCSSP(J) ,0.0,1) 

c 
c 

CHAFF  RCS 
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cal 1  symbol (4.0,0.3,0.16, CRCS ,0.0,9) 

call  number (6. 0,0. 3,0. 16, RCSCP (J) ,0.0,1) 
RETURN 
END 
C* *************************************************************** 

SUBROUTINE  SHPCLL(J) 
C    THIS  SUBROUTINE  CALCULATES  THE  POTION  OF  THE  SHIP  RCS  IN  THE 
CELL 

COMMON/MI SS/RG , THETA , AG , PI 

COMMON/ABCD/BRGM(60) ,R(60) ,SC,VS (60) ,CSE(60) ,BRGMT(60) , 
*BRGTM(60) ,RMT(60) 
COMMON/SHCL/AGS , RGS , BRGMS ( 60) , SAA, PARTS , SAAN (60) 
COMMON/SHP/SL, SB 
C 

IF  (RGS.GT. (RG-SL/( 2. 0*SC)))  THEN 
IF  (BRGMS (J) .LE. 0.0)  THEN 

BR  =  90.0  -  (BRGMS (J)  *  180. /PI) 
ELSE 

BR  =  450.0  -  (BRGMS (J)  *  180. /PI) 
END  IF 

CALL  ASPECT(J,BR,SAAR) 
S2R  =  ABS(RGS-RG) 
IF  (SAAR.EQ.90.0)  THEN 

SIR  =  0.0 
ELSE 

SIR  =  S2R  /  ABS(COS(SAAR*PI/180.)) 
END  IF 
END  IF 

IF  (AGS.GT. (AG-SL/(2.0*SC) ) )  THEN 
IF  (BRGMS (J) .GT.BRGMT(J) )  THEN 
IF  (BRGMS (J) .LE. 0.0)  THEN 

BA  =  180.0  -  (BRGMS(J)  *  180. /PI) 
ELSE 

BA  =  540.0  -  (BRGMS (J)  *  180. /PI) 
END  IF 
ELSE 

IF  (BRGMS (J) .LE. 0.0)  THEN 

BA  =  0.0  -  (BRGMS (J)  *  180. /PI) 
ELSE 

BA  =  360.0  -  (BRGMS (J)  *  180. /PI) 
END  IF 
END  IF 

CALL  ASPECT (J, BA,SAAA) 
S2A  =  ABS (AGS-AG) 
IF  (SAAA.EQ.90.0)  THEN 

S1A  =  0.0 
ELSE 

S1A  =  S2A  /  ABS (COS (S AAA* PI/ 180.) ) 
END  IF 
END  IF 
IF( (RGS.GT. (RG-SL/(SC*2.0) ) ) .AND. (AGS.GT. (AG-SL/ (SC*2 . 0) ) ) ) 


THEN 
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IF  ((RGS.LT.RG) .AND. (AGS.LT. AG) )  THEN 
IF  (S1A.LT.S1R)  THEN 
SI  =  S1A 
SAA  =  90  -  SAAA 
ELSE 

SI  =  SIR 
SAA  =  SAAR 
END  IF 
ELSE  IF  ((RGS.GT.RG) .AND. (AGS.GT.AG))  THEN 
IF  (S1A.LT.S1R)  THEN 
SI  =  SIR 
SAA  -  SAAR 
ELSE 

SI  =  S1A 
SAA  =  90.  -  SAAA 
END  IF 
ELSE  IF  ((RGS.GT.RG) .AND. (AGS. LE. AG))  THEN 
SI  =  SIR 
SAA  =  SAAR 
ELSE 

SI  =  S1A 

SAA  =  90.0  -  SAAA 
END  IF 


ELSEIF( (RGS.GT. (RG-SL/ (SC*2 . 0) ) ) .AND. (AGS.LE. (AG-SL/ (SC*2 . 0) ) ) ) 
•THEN 

SI  =  SIR 

SAA  =  SAAR 
ELSE 

SI  =  S1A 

SAA  =  90.  -  SAAA 
END  IF 

SAREA  =  SL  *  SB 
SAAN(J)  =  SAA 
S1SC  =  S1*SC 
IF  ( (RGS.GT.RG) .OR. (AGS.GT.AG) )  THEN 

PARTS  =  (SAREA  -  .5*SAREA  -  SB*S1SC) /SAREA 

IF  (S1SC.GT. (SL/2.0) )  THEN 
PARTS  =0.0 

END  IF 
ELSE  IF  (( (RGS.EQ.RG) .AND. (AGS.LE. AG) ) .OR. 
*( (RGS.LT.RG) .AND. (AGS.EQ.AG) ) )  THEN 

PARTS  =  0.5*SAREA 
ELSE 

PARTS  =  (0.5*SAREA  +  SB*S1SC) /SAREA 

IF  (S1SC.GT. (SL/2.0) )  THEN 
PARTS  =1.0 

END  IF 
END  IF 
RETURN 
END 
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C* *************************************************************** 

SUBROUTINE  ASPECT  (J, B,SAA) 
C   THIS  SUBROUTINE  CALCULATES  THE  ANGLE  BETWEEN  THE  SHIP'S  HEADING 
C   AND  THE  GATE  BOUNDARY  IT  EXITS.  IT  IS  USED  IN  SUBROUTINE  SHPCLL 
COMMON/ABCD/BRGM(60) ,R(60) ,SC, VS (60) ,CSE(60) ,BRGMT(60) , 
*BRGTM(60) /RMT(60) 
DETERMINE  SHIP  ASPECT  ANGLE  (SAA) 


C 
C 


QUADRANT  II 

IF  ( (90. 0. LE. B) .AND. (B.LT. 180.0))  THEN 

IF  ((O.O.LE.CSE(J)) .AND. (CSE(J) .LT. 90.0))  THEN 
IF  (B.GE.CSE(J))  THEN 

SAA  =  B  -  CSE(J) 
ELSE 

SAA  =  CSE(J)  +  180.0  -  B 
END  IF 
END  IF 

IF  ( (90.0. LE.CSE(J) ) .AND. (CSE(J) .LT. 180.0) )  THEN 
IF  (B.GE.CSE(J) )  THEN 

SAA  =  B  -  CSE(J) 
ELSE 

SAA  =  CSE(J)  -  B 
END  IF 
END  IF 

IF  ( (180.0. LE.CSE(J) ) .AND. (CSE(J) .LT. 270.0) )  THEN 
IF  ( (CSE(J)-B) .LE.90.0)  THEN 

SAA  =  CSE(J)  -  B 
ELSE 

SAA  =  B  -  (CSE(J)  -  180.0) 
END  IF 
END  IF 

IF  ( (270.0. LE.CSE(J) ) .AND. (CSE(J) .LT. 360.0) )  THEN 
IF  ( (CSE(J)-B) .LE. 180.0)  THEN 

SAA  =  B  -  (CSE(J)  -  180.0) 
ELSE 

SAA  =  (CSE(J)  -  180.0)  -  B 
END  IF 
END  IF 
QUADRANT  I 


ELSE  IF  ( (180.0. LE.B) .AND. (B.LT. 270.0) )  THEN 

IF  ( (0.0. LE.CSE(J) ) .AND. (CSE(J) .LT. 90.0) )  THEN 
IF  (B.GE. (CSE(J)+180) )  THEN 

SAA  =  B  -  CSE(J) 
ELSE 

SAA  =  CSE(J)  +  180.0  -  B 
END  IF 
END  IF 

IF  ( (90.0. LE.CSE(J) ) .AND. (CSE(J) .LT. 180.0) )  THEN 
IF  ( (B-CSE(J) ) .LE.90.0)  THEN 
SAA  =  B  -  CSE(J) 
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ELSE 

SAA  =  CSE(J)  +  180.0  -  B 
END  IF 
END  IF 

IF  ( (180. O.LE.CSE(J)) .AND. (CSE(J) .LT. 270.0) )  THEN 
IF  (CSE(J) .GE.B)  THEN 

SAA  =  CSE(J)  -  B 
ELSE 

SAA  =  B  -  CSE(J) 
END  IF 
END  IF 

IF  ( (270. O.LE.CSE(J)) .AND. (CSE(J) .LT. 360.0))  THEN 
IF  ((CSE(J)-B) .LE.90.0)  THEN 

SAA  =  CSE(J)  -  B 
ELSE 

SAA  =  B  -  (CSE(J)  -  180.0) 
END  IF 
END  IF 
C    QUADRANT  IV 

C    

ELSE  IF  ( (270.0. LE.B) .AND. (B.LT. 360.0) )  THEN 

IF  ( (0.0. LE.CSE(J) ) .AND. (CSE(J) .LT. 90.0) )  THEN 
IF  ( (B-(CSE(J)+180) ) .LE.90.0)  THEN 

SAA  =  B  -  (CSE(J)+180.) 
ELSE 

SAA  =  CSE(J)  -  B  +  360.0 
END  IF 
END  IF 

IF  ( (90.0. LE.CSE(J) )  .AND.  (CSE(J)  .LT. 180.0) )  THEN 
IF  (B.GE. (CSE(J)+180.0) )  THEN 

SAA  =  B  -  (CSE(J)+180.0) 
ELSE 

SAA  =  CSE(J)  +  180.0  -  B 
END  IF 
END  IF 

IF  ( (180.0. LE.CSE(J) ) .AND. (CSE(J) .LT. 270.0) )  THEN 
IF  ( (CSE(J)+180.0-B) .GE.90.0)  THEN 

SAA  =  B  -  CSE(J) 
ELSE 

SAA  =  CSE(J)  +  180.0  -  B 
END  IF 
END  IF 

IF  ( (270.0.  LE.CSE(J) )  .AND.  (CSE(J)  .LT. 360.0) )  THEN 
IF  (B.GE.CSE(J) )  THEN 

SAA  =  B  -  CSE(J) 
ELSE 

SAA  =  CSE(J)  -  B 
END  IF 
END  IF 
C    QUADRANT  III 
C    

ELSE  IF  ( (360.0. LE.B) .AND. (B.LT. 450.0) )  THEN 
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IF  ((O.O.LE.CSE(J)) .AND. (CSE(J) .LT. 90.0))  THEN 
IF  ((B-360.0) .LE.CSE(J))  THEN 

SAA  =  CSE(J)  -  (B  -  360.0) 
ELSE 

SAA  =  (B  -  360.0)  -  CSE(J) 
END  IF 
END  IF 

IF  ( (90. 0. LE.CSE(J)) .AND. (CSE(J) .LT. 180.0))  THEN 
IF  ((CSE(J)-(B-360.0)) .LE.90.0)  THEN 

SAA  =  CSE(J)  -  (B-360.0) 
ELSE 

SAA  =  (B  -  360.0)  -  (CSE(J)  -  180.0) 
END  IF 
END  IF 

IF  ( (180. O.LE.CSE(J)) .AND. (CSE(J) .LT. 270.0))  THEN 
IF  ((CSE(J)-180.0) .GE. (B-360.0))  THEN 
SAA  =  (CSE(J)  -  180.)  -  (B  -  360.0) 
ELSE 

SAA  =  (B  -  360.)  -  (CSE(J)  -  180.0) 
END  IF 
END  IF 

IF  ( (270.0. LE.CSE(J) ) .AND. (CSE(J) .LT. 360.0) )  THEN 
IF  ( (B-CSE(J) ) .LE.90.0)  THEN 

SAA  =  B  -  CSE(J) 
ELSE 

SAA  =  CSE(J)  -  (B  -  180.0) 
END  IF 
END  IF 
ELSE 
END  IF 
RETURN 
END 
C  A************************************************************** 

SUBROUTINE  CHFCLL 
C    THIS  SUBROUTINE  CALCULATES  THE  CHAFF  RCS  IN  THE  CELL 
COMMON/MISS/RG , THETA , AG , PI 
COMMON/ CFCL/PARTC , CRAD , RGC , AGC 
COMMON/RC/RCSS , RCSSP ( 60 ) , RCSC, RCSCP (60) 
CAREA  -  PI  *  CRAD**2.0 
R1R  =  ABS(RGC  -  RG) 
R1A  =  ABS(AGC  -  AG) 
IF  ( (R1R.LT. CRAD) .AND. (R1R.GT. 0.0) )  THEN 

A2R  =  ACOS (R1R/CRAD) 

R2R  =  CRAD*SIN(A2R) 

SECTR  =  A2R*CRAD**2.0 

TRIR  =  R1R*R2R 
END  IF 
IF  ( (R1A.LT. CRAD) .AND. (R1A.GT. 0.0) )  THEN 

A2A  =  ACOS (Rl A/ CRAD) 

R2A  =  CRAD*SIN(A2A) 

SECTA  =  A2A*CRAD**2.0 

TRIA  =  R1A*R2A 
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END  IF 

IF  (((R1R.LT.CRAD) .AND. (R1R.GT. 0. 0) ) .AND. 
*( (R1A.LT.CRAD) .AND. (R1A.GT.0.0) ) )  THEN 


A2A  = 

ACOS (R1A/CRAD) 

A2R  = 

AC0S(R1R/CRAD) 

R2A  = 

CRAD*SIN(A2A) 

R2R  = 

CRAD*SIN(A2R) 

SECTA 

=  A2A*CRAD**2.0 

SECTR 

=  A2R*CRAD**2.0 

TRIA  = 

=  R1A*R2A 

TRIR  = 

=  R1R*R2R 

END  IF 

IF  ((R1R.EQ.0.0) .OR. (R1A.EQ. 0. 0) )  THEN 
PARTC  =0.5 
IF  ((R1R.EQ.0.0) .AND. (R1A.EQ. 0. 0) )  THEN 

PARTC  =0.25 
END  IF 
ELSE  IF  ((R1R.LT.CRAD) .AND. (R1A.GT.CRAD) )  THEN 
IF  (RGC.LT.RG)  THEN 

PARTC  =  (CAREA-SECTR+TRIR)/CAREA 
ELSE 

PARTC  =  (SECTR-TRIR)/CAREA 
END  IF 
ELSE  IF  ( (R1R.GT.CRAD) .AND. (R1A.LT.CRAD) )  THEN 
IF  (AGC.LT.AG)  THEN 

PARTC  =  (CAREA-SECTA+TRIA)/CAREA 
ELSE 

PARTC  =  (SECTA-TRIA)/CAREA 
END  IF 
ELSE  IF  ( ( (R1R.LT.CRAD) .AND. (R1R.GT.0.0) ) .AND. 
*( (R1A.LT.CRAD) .AND. (R1A.GT.0.0) ) )  THEN 
IF  ( (RGC.LT.RG) .AND. (AGC.LT.AG) )  THEN 
IF  (R1R.LT.R1A)  THEN 

PARTC  =  (CAREA-SECTR+TRIR)/CAREA 
ELSE 

PARTC  =  (CAREA-SECTA+TRIA)/CAREA 
END  IF 
ELSE  IF  ( (RGC.LT.RG) .AND. (AGC.GT. AG) )  THEN 

PARTC  =  (SECTA-TRIA)/CAREA 
ELSE  IF  ( (RGC.GT.RG) .AND. (AGC.LT.AG) )  THEN 

PARTC  =  (SECTR-TRIR)/CAREA 
ELSE 

IF  (R1R.LT.R1A)  THEN 

PARTC  =  0.5*(SECTA-TRIA)/CAREA 
ELSE 

PARTC  =  0.5*(SECTR-TRIR)/CAREA 
END  IF 
END  IF 
ELSE 

PARTC  =0.0 
END  IF 
IF  ( (RGC.GT. (RG+CRAD) ) .OR. (AGC.GT. (AG+CRAD) ) )  THEN 
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PARTC  =0.0 
END  IF 
RETURN 
END 
C* *************************************************************** 

SUBROUTINE  SCREN3 (IFLAG3) 
C    THIS  SUBROUTINE  ENABLES  ALTERATION  OF  THE  INPUT  VARIABLES 
COMMON/ABCD/BRGM(60) ,R(60) ,SC,VS (60) ,CSE(60) , BRGMT(60) , 
*BRGTM(60) ,RMT(60) 
COMMON/WD/GW , VW 
COMMON/MISS/RG , THETA , AG , PI 
COMMON/RC/RCSS,RCSSP(60) ,RCSC,RCSCP(60) 
INTEGER  I 


IFLAG3=0 
100  WRITE(*,*) 
WRITE (* 
WRITE (* 
WRITE (* 
WRITE (* 
WRITE (* 
WRITE (* 
WRITE (* 
WRITE (* 
WRITE (* 
WRITE (* 
WRITE (* 
WRITE (* 
WRITE (* 
READ  (* 
IF 


*) 
*) 
*) 
*) 
*) 
*) 
*) 
*) 
*) 
*) 
*) 
*) 
*) 
*) 


(1) 
(2) 
(3) 
(4) 
(5) 
(6) 
(7) 
(8) 

(9) 
(10) 

(11) 
(12) 

(13) 


' INPUT 
I 
(I.EQ.l)  THEN 

WRITE (*,*) 

READ  (*,*) 

IFLAG3  =  2 

GO  TO  100 
ENDIF 
IF  (I.EQ.2)  THEN 

WRITE (*,*)  'INPUT 

READ  (*,*)  RANGE 

IFLAG3  =  2 

GO  TO  100 
ENDIF 
IF  (I.EQ.3)  THEN 

WRITE (*,*)  'INPUT 

READ  (*,*)  CSE(l) 

IFLAG3  =  2 

GO  TO  100 
ENDIF 
IF  (I.EQ.4)  THEN 

WRITE (*,*)  'INPUT 

READ  (*,*)  SC 

IFLAG3  =  2 

GO  TO  100 


MISSILE  BEARING' 
MISSILE  RANGE' 
SHIP"S  HEAD' 
SCREEN  SCALE' 
MISSILE  RADAR  BEAMWIDTH' 
RANGE  GATE' 

MAXIMUM  RADAR  CROSS  SECTION  OF 
MINIMUM  RADAR  CROSS  SECTION  OF 
RADAR  CROSS  SECTION  OF  CHAFF' 
TRUE  WIND  DIRECTION' 
TRUE  WIND  SPEED' 
ALL  PARAMETERS' 
RETURN  TO  MAIN  PROGRAM' 
SELECTION' 


SHIP' 
SHIP' 


'INPUT  MISSILE 
BRGM(l) 


BEARING 


INITIAL  MISSILE  RANGE 


SHIP"S  HEAD' 


SCREEN  SCALE' 
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ENDIF 

IF  (I.EQ.5)  THEN 

WRITE (*,*)  'INPUT  MISSILE  RADAR  BEAMWIDTH' 

READ  (*,*)  THETA 

IFLAG3  =  2 

GO  TO  100 
ENDIF 
IF  (I.EQ.6)  THEN 

WRITE (*,*)  'INPUT  RANGE  GATE' 

READ  (*,*)RGATE 

RG=RG/SC 

IFIAG3  -  2 

GO  TO  100 
ENDIF 
IF  (I.EQ.7)  THEN 

WRITE (*,*)  'INPUT  MAXIMUM  RADAR  CROSS  SECTION  OF  SHIP 

READ  (*,*)  RCSSMAX 

IFLAG3  =  2 

GO  TO  100 
ENDIF 
IF  (I.EQ.8)  THEN 

WRITE (*,*) 'INPUT  MINIMUM  RADAR  CROSS  SECTION  OF  SHIP' 

READ  (*,*)  RCSSMIN 

IFLAG3  =  2 

GO  TO  100 
ENDIF 
IF  (I.EQ.9)  THEN 

WRITE (*,*)  'INPUT  RADAR  CROSS  SECTIONOF  CHAFFf 

READ  (*,*)  RCSSC 

IFLAG3  =  2 

GO  TO  100 
ENDIF 
IF  (I.EQ.10)  THEN 

WRITE (*,*)  'INPUT  TRUE  WIND  DIRECTION' 

READ  (*,*)  GW 

IFLAG3  =  2 

GO  TO  100 
ENDIF 
IF  (I.EQ.ll)  THEN 

WRITE (*,*)  'INPUT  TRUE  WIND  SPEED' 

READ  (*,*)  VW 

IFLAG3  =  2 

GO  TO  100 
ENDIF 

IF  (I.EQ.12)  IFLAG3=1 
IF  ((I.LT.l) .OR. (I.GT.13) )  THEN 

WRITE (*,*)  'INPUT  AGAIN' 

GO  TO  100 
ENDIF 
RETURN 
END 
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Deployment^!  chaff  in 
centroid  mode  against 
anti-ship  missiles  using 
a  variable  azimuth  and 
elevation  launcher. 
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